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SYITOPSIS 


A Thesis entitled "Some Contributions 
to Mathematical Programming” submitted in 
Partial Fulfilment of the Requirements 
for the Degree of Doctor' of Philosophy 
by Purshotam Lai Jolly to the Department 
of Mathematics, Indian Institute of Technology 
Kanpur « 


This thesis discusses pseudo-monotonic programming 
and its applications, the sensitivity analysis for linear 
complimentarity problem with applications, boxnded variable 
linear programs and the properties of assignment polytope 
and traveling salesman polytope* 

Chapter 1 is introductory in nature and is intended 
to bring out the contributions made in this thesis in their 
proper perspective* 

In Chapter 2, a linearization technique to solve 
programming problems with pseudo-montonic objective 
functions is given* The algorithm solves the pseudo- 
montonic program by solving a sequence of finite number 
of linear programs over convex or nonr-convex sets (e*g« 
set of discrete points). The applicability of the 
algorithm is demonstrated by considering pseudo-monotonic 
integer programs (all integer and mixed integer case), 
zero-one programs, transportation, assignment, traveling 
salesman problems, set covering and set. partitioning 
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problems and programming problems with linear 
objective f-unctions. 

In Chapter 3, we discuss the sensitivity analysis 
of linear compliment arity problem (LCP) 

w - Mz = q 

, / V / 7 •) M is nxn 

where w= Cwj^jWg, . , .jW^^; , z = (,z^,Zg, 

matri^c and q is n-vectoi’? using Lemke>s 

pivot algorithm [^"Bimatrlx Equilibrium Points and 

mmmrn 

Matheraatical Programming”, Man.Soi.jlljlSCbjpp^SSVSS * > 

discussing the changes in vector q, changes in one 
more columns of matrix M and simultaneous changes in 
vector q and matrix M* Sensitivity analysis of oon 
quadratic programming problem is discussed through 

LCP. As an application of the procedure, a simple 

. 1 fractional 

computational scheme is devised to solve lineax 

functional programs and quadratic fractional functional 
programs. 

In Chapter 4, we propose two methods to solve 

bounded variable linear programs. The approach use 

in both these methods is essentially a relaxation app 

In Chapter 5, we study some aspects of assiguiaent 

polytope (APj^) by using an elementary approach, ^ 

* 0 ^ 

particular, we study the adjacency of the vertices 

* *1 ]|* C* jj - Tld s 

assignment polytope, the structure of some specia 

of its faces, the diameter of AP^ and the edge- connectivity 
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of the graph of AP^* 

In Chapter 6^ w© make an attempt to study the 
adjacency of the vertices of the travelijag salesman 
polytope, A necessary condition for non-adjacency 
of t-wo vertices of traveling salesman polytope is given* 
A sufficient condition for non-adjacency is also 
established. However, the necessary condition is not 
sufficient and the sufficient condition is not necessary, 
X'le also formulate some adjacency and non-adjacency rules 
to generate adjacent and non- adjacent vertices to a 
given vertex on traveling salesman polytope and discuss 
their applications. 
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INTRODUCTION 

There are many problems in government, defence 
and industrial organisations which are concerned with 
picking up of a best alternative from amongest a large 
number of available alternatives. One of the important 
class of such problems is the allocation of scarce 
resources in an optimum manner. Some such problems 
are to plan the economy of a nation in a most efficient 
manner, deployment of the armed forces in a way so as 
to maximize the country* s chance of winning a war, to 
maximize the profit of some organisation or to minimize 
the waste of some products etc. These problems are 
so important that no organisation can afford to ignore 
them. The formulation of such problems and their methods 
of solutions gave rise to a new class of optimization 
problems which cannot be solved by the classical methods 
of optimization based on differential calculus or 
calculus of variations. Among them one is the class of 
Mathematical Programming Problems , This thesis makes 
some contributions to this class of problems. 

The main purpose of this chapter is to give a 
brief survey of theoretical and computational developments 
in mathematical programming which are relevant to this thesis* 
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The purpose of this survey is to hring out the 
contributioB more clearly, rather than being exhaustive. 

To make the survey more precis© and understandabXe the 
present chapter is divided into five sections. Section J.*!. 
gives preliminaries and definitions. Section 1.2 gives 
the classifications of mathematical programs. Section 1,3 
gives their brief history of development. Section 1,4 
discusses relevant computational techniques used in this 
thesis and in section 1.5, we discuss contributions 
of the thesis in brief, 

1,1 Pr eliminaries and Definit^or^ s 

For the sake of completeness, we give here 
some of the standard definitions. 

Let f ; X R where X C R^ and x^,x^ be 
any two points of X, Let x° £ X be such that 

= (l-.oc)x^ + oCx2 for 0 < °( < 1, Clearly x° G X if 
X is a convex set. Then we say f is [ 112] 

(i) strictly convex : if X is a convex set and 

f(x°) < Cl-'=c)f(x^) + o(f (x^) ^ 

(ii) convex : if X is a convex set and 

f (x°) < (1-oC) f(x^) + oCf(x^) 

(iii) Dseudo~ convex : if f is continuously differentiable 
on an open set containing X and 

Vf(x^) (x^-x^) > 0 fCx^) > f<x^) 

mm 
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Civ) explicit quasi'»» conve x s if Z is a convex set and 

f(x°) < max ij^(x^), f(x^)Jl» , f(x^) £(x^)t ‘=< 0 I 

(v) quasi- convex : if Z is a convex set and 

f(x°) < max (x^),f (x^)^ 5 or 
equivalently, f is quasi- convex if it is 
continuously differentiable on an open set 
containing X and 

f(x^) < f(x^) Vf (x^)(x^-x^) < 0 

mm ^ 

Likewise we say that f is strictly concave 
(cancave, pseudo-cancave, explicit quasi- concave, 
quasi-concave) if -f is strictly convex (convex, 
pseudo-convex, explicit quasi-convex, quasi-convex) 
on Z* 

(vi) paeudo-monotonic : if f is both pseudo-convex 
and pseudo-concave. 

(vii) quas i- montonic ; if f is both quasi-convex and 
quasi- concave. 

Below we give some definitions concerning 

convex polytopes [41,71]. 

(i) convex polyhedron : By a convex polyhedron we mean 
solution set Z of the following system of linear 
constraints 

Ax = b, X > 0 

where A,b are given mxn and mxl matrices respectively 
and X is a nxl column vector. 
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(ii) convex polvtope : A bounded convex polyhedron 

is called a convex polytope* 

Ciii) Extreme nnint (vertex) t A point X 0 X is called 
an extreme point (vertex) of X iff it cannot 
be expressed as a convex combination of any 
two distinct points of X* 

A convex polytope can be regarded as the 
convex hull of its finite number of extreme points. 

(iv) adj acent extreme points t Two extreme points 
x^ and x^ S X are called adjacent (neighbouring) 
extreme points on X (or adjacent simply) iff 
every point on the line segment joining x and 

O 

X is uniquely expressed as a convex combination of 
extreme points of Xm 

(v) edge t The line segment joining any pair of adjacent 

extreme points of X is called an edge of X. 

1 2 

Thus X ,x are adjacent on X iff the line segment 
1 2 

joining x and x is an edge of X, 

We adopt the following notation to indicate 

T 2 

the line segment joining x"^ and x • 

[x^,x^] = •^x : X = °<x^ + (l<-o()x^, 0 < <=< 

(vi) supporting hvperplane ; Consider a boimdary 

point X?- of convex polytope (or convex set)X. 

Then cx = z where c = (cn,c , .*.,0 ) is called a 

1 1 

supporting hyperplane at x-^ if cx = z and^^hole 
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Of X lies in one closed half space produced by 
the hyper plane i.e, cu> z # uEXor 
cu < z ¥ u 6 X, 

(vii) face : Face of a polytope is the intersection 

of the polytope with some of its supporting hyperplane, 
(viii) permutation matrix i A nxn zero~one matrix 

is called a permutation matrix if each of its rows 
and columns contains exactly one 1. 

The above definitions are illustrated in the 
following figures. 



We now giveafOw definitions of graph theory [ V8]« We 
assume the definitions of a graph, its nodes and of 
edges. 

(i) valency of a node : 7%. valency of a node is the 
number of edges incident to it., 
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(ii) connected graph : A graph G is ealled connected 
if each node of it can he reached through edges 
and nodes of it* 

(iii) p, ath : A path in a connected graph G is a connected 

subgraph of the se-me having two nodes of valency one 
and all other nodes of valency two* 

(iv) length of a path : The length of a path is the 
number of edges in it. 

(v) circuit t A circuit in a graph is a closed path 
having all nodes of valency two, 

1,2 Optimization Problems 

The general optimization problem can be 
described as follows : 

Let A be an arbitrary set and B a set with 
a transitive relation T defined in it. Let f : A B, 

Then the optimization problem is to find an element 
x° $ A, if it exists, such that 

(i) f (x) T f (x®) ¥ X e A - <.<=> 

or (ii) f(x) •¥ x G A - 

Usually, the set A consists of all possible 
decisions or policies under consideration, the set B is a 
subset of real numbers, R, and f is a criterion or ob.iect:^ve 
function i,e, a rule that associates to each decision in A 
a unique element in B, 
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With any optiinization problem ’we are primarily 
concerned with the following types of questions. 

(i) Does such an exists ? 

(ii) If such an x° exists, is it unique ? 

Ciii) If such an x° exists, is there a 

procedure to find it out ? 

Optimization problems are subdivided into classes 
of problems either based upon the model or upon the 
solution procedure. 

In the following, we give the various problems 
which today constitute Mathematica l Programming Problems . 

1,2*1 Mathematical Programs 

A mathematical program is a particular case 

of the optimization problem. Here set B stands for the 

set of real numbers, T is the usual ordering C < , 2 ^ 

in the set of real numbers and A known as^constraint set, 

n 

is a subset of B which is usually of the form : 

A =-^ e R^l g^(x) ( 

where R, i = l,2,..,,m, are called constraint 

functions or simply constraints . 

Here below we classify the mathematical programs. 
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1*2.2 |::;i,asslfication of Mathematical grograms. 

Broadly speaking} the mathematical programs can he 
classified into two main categories, namely 

(A) Stochastic Math ematical Pr ograms : The programs 

in which the objective function f and/or the 
constraints i = 1, involve elements 

of chance are called stochastic mathematical 
programs* 

(B) Deterministic Mathematical Programs s The programs 
which are not stochastic in nature are called 
deterministic mathematical programs* 

We, in this thesis, concentrate on deterministic 
mathematical programs only, 

1,2.3 Classification of Deterministic Mathematical Programs 

Depending upon the nature of f and gj_’s ( i = l}*,m) 
deterministic mathematical programs can now he broadly 
classified as under* 

B1 Linear Programs : Mathematical programs in which 
the objective function f and the constraints 
g^s(i = l,,,,,m) are linear are called linear 
programs* Another feature of these programs is 
the imposition of non-negativity restrictions 
over the variables (j = 1,2, ...jn). Linear 
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programs have large number of applications and 
are widely studied programs. They can be written 
as 

max f(x) f cx 
subject to 

Ax = b 
X > 0 

where c is an n-vector, b an m-vector and A is 
an mxn matrix. 

Linear programs have been studied in more general 
forms viz, 

Bll Sem i- I nfinite Linear Programs : The linear 

programs in which m becomes infinite (Gharnes, 

Cooper and Kortanek [32,33,34,35] ). 

B12 Infinite Linear Programs t The linear programs 
in which both n and m become infinite (Duffin 
and Kartovitz [ 47 ] ) , 

B2 N on-Linear Programs : A mathematical program 

which is not linear is called non-linear program. 

Thus non-linearity can be introduced by taking 
f and/or one or more of g^*s (i = l,...,m) as 
non-linear functions. There is a considerable interest 
in the study of non-linear programs because of their 
challenging nature. Some of them have been so 
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extensively studied, that they have acquired the 
status of a separate subject of their o\'m, A 
few important classes of these programs are 
listed below : 

B21 Convex (Concave )Programs ; In these programs we 
minimize (maximize) convex (concave) - like 
functions subject to convex/concave - like 
constraints* As a linear function is convex 
as well as concave, it can be treated aso-convex 
as well asaconcave program, Quadratifi programs 
in which the objective function f is quadratic 
i,e, f(x) = 3 C + cx and (i = l,,.,,m) are 

linear is another example of such programs. The 
variables x^ (j = 1, ,,,,n) are assumed to be 
non-negative. Depending upon the natm’e of 
Q, a quadratic program becomes convex or concave. 
The other programs in this class include 
piece-wise convex programs (Zangwill [158] , 
Charnes, Cooper and Kortanek [35] ) , pseudo- 

convex programs (Mangasaraian [ill] ), generalised 
pseudo-convex programs (Gupta and Bhatt [73] ), 

explicit quasi-convex programs (Bela Martos [H^] )j 
quasi-convex programs (Arrow and Enthoven [5] ). 
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3P2 Fractional Programs i In these programs, f ^ q 
where p ; ^ E and q ; R^ R. If p,q and 

(i = 1, are all linear functions 

then the program is called linear fractional 
program otherwise it is the case of non-linear 
fractional program. 

Few important papers concerning fractional 
programs are (Charnes and Cooper [ 3l] , Dorn [45} , 

Ivanti Swarup [ 141, 142,144,145, 146] , Aggarwal [1,2], 
Dinkelhach [44] , Sector [15,16] , Gogia [65] 

B23 Indefinite Programs : If f is the product of two 
functions p and q, the mathematical program is 
called an indefinite program ( Sector [16] , Kanti 
Swarup [ 145,146] etc, ), 

B24 Integer Programs s A mathematical program in 

which the variables are restricted to be integers 
is called an integer program. If all the variables 
are restricted to be integers, it is called an all 
integer program otherwise a mixed integer program 
(Gomory [66,67]). If the objective function f and 
constraint set (i = 1, ,,,,m) are all linear, 

"‘fcCk 

it is sometimes referred ^as«.linear integer program. 
Likewise we may have fractional integer programs 
and quadratic integer programs. 
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B25 Zer o-One Pro£r am_s t A mathematical program 
in which variables are restricted to he zero 
or one is called zero-one program (Balas [8] , 

Geoffrion ^60] | Glover [ 6l] ), 

B26 Geometric Programs : A mathematical program in 

which 

o ri a„. . 

fCx) = 2 c . rr ^ 

j=l i=l ^ 

and the constraint set A is given hy 

A e R^jg^(x) < 1, k = l,2,,,,,m^ 

^ Si • • 

where g, (x) = 2 c, . TT x. = 1 

k 1 > 

where represents the number of terms in the 
objective function, and denotes the number of 
terms in the kth constraint, c^^’s, c^^’s and 
akij*s are real numbers (Duffin [46] , Duff in and 

Peterson [48] , Passy and i/Jilde [122] ), 

1*2 .4 Programming Problems in Abstract Spaces : 

It is not always that underlying set A is a 
subset of finite dimensional Euclidean space R^, A 
mathematical program in which underlying set A is not a subset 
of finite dimensional Euclidean space is included in a class 
of programming problems in abstract spaces# Such programming 
problems have been studied in locally convex space 
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(Hurwicz [82] , Hanson [77] ), Hilbert space (Kelley and 
Thompson [94], ), pre Hilbert space (Blura [26] ), 

The following problem in real Banach space has been 
extensively studied (Varaiya ,[ 153]), 

max f(x) 

subject to g(x) £ Ay, x £ 

where A ^ d X and 
g jX-^Yandf ;X-*R, 

For further results, refer to (Ritter [131] , 

Grignard [70] , Neustadt [ 120] j Pchenichney [123] ), 

1,3 B rief H istorical Review 

Optimization problems have been of keen interest 
to mathematicians, engineers and physical scientists 
since long. Since the end of word war II there has been 
a rapidly increasing interest in business, management, 
defence and scientific world in the optimization problems. 
Although Kantorovitch reported some results in 1939, the 
real impetus came only in 1947 when G,B,Dantzig gave a very 
powerful technique called Simplex Method for solving 
linear programming problems. Although the simplex method was 
developed in 1947, it was generally not available till 1951, 
when the cowless commission monograph no, 13 was published. 


Ay dX, X and Y are real Banach spaces, 
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Soon after that attention was paid to linear programming 
problems with special structoe like transportation, 
assignment and net work problems etc, and were solved 
successfully. Charnes and Cooper [30] , Ford and 

Fulkerson [54,55] and H,W,Kuhn [ 98j made significant 
contributions in this direction. The invention of high 
speed digital computers added more incentive to the 
developing interest in linear programming problems* David 
Gale, H*l'J*Kuhn and A,\f .Tucker made important contributions 
in developing the duality theory for linear programming* 
C.E.Lemke made a major contribution in this direction 
by developing a dual simplex technique [ 103 ] for 
solving linear programs. 

There are many practical problems which cannot 
be represented by linear programming models* Hence the 
need for more general mathematical progra^mming models 
led to the growing interest in non~ linear programming 
problems. However, the work in this direction started 
simultaneously. The first set of necessary and sufficient 
conditions for the existence of optimal solutions to nonr- 
linear programming problems came as early as 1951 when 
H.W.Kuhn and A,¥*Tucker [99] published an important paper 
’’Non-linear Programming”, Their work laid the foundation 
of a great deal of later work in non-linear programming* 
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However, before Kuhn and Tucker, Fritz John [ 88 ] 
derived a set of necessary conditions for the case of 
inequality constraints alone. The other generalizations 
include the work by Arrow, Hurvicz and Uzawa L 6,7 ] , 
and Kangasarian [111,113] ate. 

The work on integer solutions to linear programming 
problems also began early in the development of the subject 
of mathematical programming. Dantzig, Fulkerson and 
Jolinson [43] published a paper in 1954 on^^raveling 
salesman problem which was concerned with the subject of 
integer solutions to linear programming problem. In 1958,'’ 
Gomory developed a cutting plane technique, a systematic 
computational method, for integer solutions to linear 
Pi’ogramming problems. In 1958, he settled ^all integer 
case [66] and in 1960 a mixed integer case [57] , 

The methods converge in finite number of iterations and 
are extensions of dual simplex technique for linear 
Programs, Later on, a number of other methods have 
appeared. For example, Ben-Israel and Gharnes [24]., 

Young [ 166 , 167 ], Fred Glover [62,63,64] , Land and Doig [lOO] 

Special attention has been paid to the zero-one 
integer programs. Work of Egon Balas [ 8 ] , Geoffrion [ 60 ] 
and Glover [ 61 ] is significant contribution in this 
direction. Their methods use implicit enumeration 
procedure which is essentially a branch and boixid procedure. 



16 


There are some zero~one programming problems 
which, because of their special structure, have attracted 
considerable attention of the research workers in the 
recent past, Hiapsack problem, the traveling salesman 
problem, set covering and set partitioning problems are 
few of the important zero-one programs, A few of the 
important papers in this direction are by BellmanlSl] , 
Bellmore and Neinhauser [22 ] , Eastman [^9] , Little, 

Murty, Sweeney and Karel [109} ,Held and Karp [ 79,80 ], 
Bellmore and Ratliff [ 23 ] , Garfinkel and Nemhauser [ 57 ] , 
Lemks, Salkin and Spielberg [ 107 ] , Balas and 
Padberg [9,10]. etc. Cutting plane techniques [ 23 ] 
and implicit enumeration procedures [lO7|3C0]are shown 
to be very effective for the solution of these problems. 

The work in quadratic programming also started 
simultaneously. The most well kno-vjn methods are by 
Beale [13,14] , Van de Panne aiid ¥hinston [150,151,152 ], 
Houthakker [ 81 ] and Wolfe [155] and thus the quadratic 
programming problems can be solved by simplex like 
procedures. The work in linear programming and quadratic 
programming has attained sufficient maturity by this 
time. Efforts have been made to provide a unified theory of 
for studying linear programs and quadratic programs. This 
link has been established by the complimentary theory of 
linear programming, Lemke>s complimentary pivot method [104 ] 
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is a major advance in this direction. Cottle, Dantzig 
and Eaves have also their major share in developing 
the theory of complimentarity. 


1.4 


In Chapters 2,3 and 4 we are concerned id-th 
some algorithms for mathematical programs. In this 
section, we give a brief survey to provide a background 
of the methods used there* 


1 • G omorv* s Cutting Plane Algorithm 

Cutting plane methods have been devised to solve 
integer linear programs. Using the optimal solution 
of the linear program (without integer restrictions), 
if not already an integer solution, a constraint, called o- 
Cut , is obtained which is such that no integer point 
is cut off and the LP solution is cut off from the 
original feasible region. The UP is again solved andAadditional 
cut is introduced at each iteration if necessary, 

Gomory [66,67] has derived cutting planes for all integer 
solutions and mixed integer solutions to linear programs. 


2, Branch and Bound Methods 

There are many combinatorial optimization 
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problems which are very difficult to solve* Some such 
problems are traveling salesman problems, scheduling 
problems, Integer programs etc* Branch and bound 
methods provides a strategy to solve moderate size 
problems of the above category* The method essentially 
consists of partitioning the set of feasible solutions 
into a finite number of subsets and then applying the 
branch and bound strategy itself to each of the so called 
candidate problems. The rules which give rise to the 
subsets are called br anc h Ing rule s * An attempt is also 

made to obtain a lower (upper) bound for each of the 
candidate problems. It is obvious that branch and bound 
methods are problem dependent in the sense that branching 
ruj,es and bounding rules depend very heavily on the 
problem Itself, The branch and bound strategy has been 
used very effectively for the following problems s 

(i) 

(ii) 

(iii) 

(iv) 

(v) 


Traveling Salesman Problem (Little, Murty, 

Sweeney and Karel [ 109 ] ) 

Scheduling [ 36 ] 

Zero-One Programs (Balas [8 ] j Geoffrion [ 60 ] 

r* 

and Glover [ 61 } ) 

Set Covering Problem (Lemke,Salkin and SpielbergjlO^■]^ 
Set Partitioning Problem (Garfinkel and 
Nemhauser C 57 ] ) 
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3 • LemkeJ_s_ A lgorithm 

LerakeJs pivot algorithm starts with an almost 

complimentary feasible solution which is obtained by 

introducing an artificial vector z = (-Ij-lj . » . «“!) 

o 

into the basis. At each stage the algorithm introduces 
a vector which is complimentary to the vector just 
dropped from the basis. All the intermediate basic 
vectors obtained in this algorithm are almost complifflentary 
feasible basic vectors. If at some stage of the 
algorithm a complementary feasible basis vector is 
obtained, it is a final basic vector and the algorithm 
terminates. For details see [119,124], 

4, Dinkelbach^s ilethod 

Dinkelbach [44] developed a technique to solve 

the following non~ linear fractional program 

min N (x )/D (x ) 

xes 

where S is compact and connected subset of and N(x) 

and D(x) are real valued continuous functions on S 

such that D(x) > 0 for all x 6 S. 

The technique depends upon the following theorem : 

Theorem 1 i q = N(x^)/D(x^) = max W(x)/D(x) iff 

° ° xes 

F(q-) = F(q ,x ) = max N(x)~q_DCx) = 0 
o ° ° xes ° 
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The algorithm is a.s follows i 

Step , 1 ; Set q = (Initially i = 0) 

Ste^J. s Solve F(q. ) = min f n(x) - q D(x)\ 

^ xes ^ ^ J 

If F(q. ) = 0 for X = then x^ solves the fractional program, 

step 3 : If F(q.) 0, compute q.,T = and 

1 1+1 

go to step 2 X'/ith i replaced hy i+1, 

5 • Bounded variahle Linear Program (BVLP) 

BVLP is of the following type 

max cx 
subject to 

Ax = h 

< u^. , j e j 

X, > 0 for all j 

3 

The simplex method is surely applicable to 
such a problem but each bounded variable will have 
to be considered as a separate constraint which results 
in an increase of the basis size. Taking into consideration 
the simple nature of the upper bound constraints f 
Dantzig [ 40 ] extended the simplex method itself to 
solve BVLP without increasing the size of the basis* The 
details of the method are available in most of the standard 
text books on linear programming a,g, [ 76,119,148 1 , 
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Contributions of the Thesis 

I'Je have givs In brief the contributions of 
the thesis* 

1 ,5 . 1 Pseudo- mo notonic Programs 

It is often found in the literature that when 

an algorithra is devised for linear objective functions, its 

possible extensions to fractional functions and other 

non-linear functions is sought. Linear fractional programs 

where the objective function is a ratio of two linear 

functions are of great interest. If the numerator 

represents the total wastage and the denominator represents 

*1 1^01 

the cost of production then ^objective ma.y be to minimize 
the wastage per unit production cost. Similarly if 
numerator represents the total cost and /^denominator 
represents the volume changes then we may be interested 
sometimes in minimizing cost per unit voluiue. Methods 
for solving linear fractional programs are given by 
Charnes and Cooper [ 31 ] , Dinkelbach [44] and Bela 
Mart os [114] , 

Linear fractional functions are special case of 
pseudo-monotonic functions methods for which, in general, 
are not available. In Chapter 2, we develop an algorithm 
for this more general class of functions. The method is a 
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linearization teclinique and provides solution to a 
pseudo-monotonic program in a finite numler of steps ty solving a 
sequence of programs with linear objective functions. 

The method is based upon a theorem which is an extension of 
the following theorem due to Kortanek and Svans [ 97 ] • 

Theorem 2 : Let f be a pseudo-concave function 

on a closed, convex set X and x' £ X, Consider the two 
programs ; 

I max f(x) , subject to x £ X. 

II max vf (x )x, subject to x £ X. 

^ * 

Then x is an optimal solution of Program I iff x is an 

optimal solution of Program II* 

¥e extend this theorem for discrete set X 

which is non- convex. 

Our algorithm for pseudo-monotonic programs 
solves all those pseud o-rnonotonic programs for which 
a solution procedure exists in case of linear objective 
functions* In particular, we have solved pseudo-monotonic 
integer programs (all integer and mixed integer), pseudo- 
monotonic zero-one programs, pseudo-monotonic transportation 
and assignment problems, pseudo-monotonic traveling 
salesman problem, pseudo-monotonic set covering and set 
partitioning problems. 
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1.5,2 Mnear _CpmDlimentarity Problem 

The Linear Coiaplimentarlty Problem (LCP ) is 
•w-Mz = q 

(LCP) w > 0, z > 0 

\/z = 0 

where w»z and q are n- vectors, M is an nxn matrix. 

There are many problems which can be formulated as 
LCP, The linear programming problems, convex quadratic 
programming problems, the problem of computing an 
equilibrium pair of strategies for bimatrix games are 
some of the problems which can be formulated as LCP, 

For example, linear programming problem can be formulated 
as LCP by combining primal and dual problems and the 
complimentary slackness property. The primal problem 
and its dual are 

min cx 
subject to 

(PP) Ax > b 

x> 0 

rp 

max by 
subject to 

,T . T 
A y < c 

y> 0 


(DP) 
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If u > Oj V > 0 are surplus and slack variables in (PP) 
and (DP) respectively then at the optimal solution 


(GS) v'^x + u^y = 0 


Taking 


r 


r 


r T ' 


T " 

V 


X 


c 


0 -A 


, z = 

1 

1! 




u 

I 

_ 

1 

i 

-b 

1 

A ■ 6 



1 

— 


the (PP)j (DP) and (CS) ctin be put as 


w-Mz = q 

It is to be noted that M is a square, asymmetric, 
positive semi-definite matrix (PSD) of order (m+n). 

Details of the problems which can be form-ulated 
as LCP are available in [ 119,124 ] . 

Dantzig and Cottle [37,38] and Lemke [104 ] 
have given methods to solve LCP of which Lemke^s method 
is applicable to a large variety of problems. Various 
conditions have been imposed over the matrix M under which 
Lemke’s method either finds a solution or shows that there 
is no solution to the problem. If the matrix M is PSD, 
which is so in case of LCPs corresponding to linear 
programs and convex quadratic programs, Lemke>s method 
finds a solution if there is any i.e. if the problem is 
feasible and unbounded, otherwise it will terminate in 
ray - showing that there is no solution. 



Details of tbs algorithm are available in [ 119,124 ] • 
Also various kinds of matrices H have been defined for 
which it has been shown that LCP has a solution [ 105 ] • 
In Chapter 3, we discuss the sensitivity 

analysis of LCP considering changes in elements of q 
and M using Lemke’s method* As an application of the 
sensitivity procedure, it is shomi that linear fractional 
functional programs and quadratic fractional functional 
programs can be solved considering sensitivity analysis 
of corresponding LCPs with changes in q only, A simple 
computational scheme is devised to solve such problems# 


1*5 #3 Bounded Variable Linear Programs 


Dantzig has proposed a method to solve BVLP 
restricting the basis size to m only. The methods has the 
following drawbacks : 

(i) Though the criteria for the selection of the 

vector to enter the basis is the same as in simplex 

method, the criteria for the vector which leaves the basis 

is decided by the row in which min 0 j-O-jU I occurs 

/ Xo. h 1 ^ 3-J f u “Xg. 

where 0 t = min W -Si- , yi i > 0 W , 0o = min J 

^ 1 i yij ' / ® 1 \-yi3 


and u. is the maximum level (upper bound) to which x. 
can be raised# This criteria is pretty involved. 
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(il) There are many problems where optimal 
solution to linear program ( without upper boutids) 
violates only a few upper bound constraints* Dantzig^s 
method does not give us any information about this* 

In Chapter 4, we propose two methods to 
solve BVLP, The proposed methods take care of above 
drawbacks, to a certain extent. One of our methods 
makes use of Dantzig’s method while the other is 
completely independent of it. Both the methods start 
with the optimal solution of linear program (without 
upper bounds). Our approach is essentially a relaxation 
approach i.e. we rela:}c some constraints at a time and 
then introduce them sequentially to arrive at the 
solution of BVLP, 

1 #5 ,4 Assignment Polvtope 

Assignment polytope XArhich arises from the 
consideration of the assignment problem is the solution 
set of the following system of linear equations 

i ~ 1, 2, • • , ,n 

0 — l,2j,,,,n 


n 

2 

X. . = 

1 

3=1 

13 


n 

2 

X. . = 

1 

i=l 

13 



^3 i 

0 
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The extreme points of the assignment polytope 
of order n are precisely the permutation matrices 

of order n [ 92 ] , It is this considera/tion which 
allows the replacement of the condition = 0 or 1 
in the assignment problem b5'- Xj . > 0 and regards AP 
as convex hull of permutation matrices of order n« 

In Chapter 5, we discuss the properties of 
the assignment polytope. Our approach is completely 
elementary. In particular, study the adjacency of 
the vertices of the assignment polytope, the structure of 
some special kinds of its faces, the diameter of AP^. 
and the edge-connectivity of the graph of AP^. By 
identifying the vertices of AP^, which are permutation 
matrices of order n,with permutations, we prove the 
following major result. 

Theorem 3 i Two permutations ir and c are adjacent 
on Ap^ iff cfTT"^ is cyclic, 

¥e have further proved that the diameter of 
the assignment polytope is tv/o and that it is edge-W(n)- 
connected where 

B(n) = 2 n (r-l)l 

r=2 

is the number of vertices adjacent to a given vertex on 

APj^, 
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1 • 5 , 5 !SLaxeJU nj^;.. Sales P olytope 

If in the above system of linear equations 
we further impose a condition that the solution set of x- .*s 

J~ J 

(such that = 0 or 1) forms a tour where by a tour 
(or tour solution) we mean a solution of the type for 
which = ... = = 1 and 

zeroes elsewhere, then the related problem is referred to 
as traveling salesman problem. The associated convex 
polytope is called the traveling salesman polytope (IPj],) 
of order n which is the convex hull of cyclic permutation 
matrices of order n, (By a cyclic permutation matrix 
we mean the permutation matrix such that the position 
of one>s give a tour solution. Cyclic permutation 
matrix is identified by a cyclic permutation of order n 
i,e, the tour on n symbols), can be regarded as 

collapsed polytope [ 118 1 of order AP^^, 

In Chapter 6, we make an attempt to study 
the adjacency of the vertices of TP^, The approach used 
for APj^ in Chapter 5 is extended in the case of TP^, A 
necessary condition for two tours to be non- adjacent 
on TP^ is given. A sufficient condition for two tours 
to be non- adjacent on TP^ is also given* Counter examples 
are given to show that neither the necessary condition is 
sufficient nor the sufficient condition is necessary. Some 
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adjacency rules and non- adjacency rules a.re formulated on 
TPn which generate tours adjacent and non-ad jacent 
respectively to a given tour. These rules can he 
used to develop some good heuristic procedures for 
solving traveling salesman problems, finally it is 
sho\-jn th8.t is hamiltonian. 



CHAPTER - 2 

PS^^QrMQ^ON IC P ROGRAMI^BTG Aro _lTS APPLI CATIONS 
2 , 1 ;[HTR pD UCTIOH 

It is well known that for a large class of 
non-linear programmes, there exists a linear programme 
such that the optimal solutions of the two prograiiimes 
are same. This is essentially an outcome of the Kuhn 
Tucker theory. However, such a linear programme is 
known only in terms of the optimal solution of the 
non-linear programme. For example j Kortanek and Evans 
proved the following theorem t 

Theor em 1 (Kortanek and Evans [97] ): Let f be a 
pseudo-concave function on a closed, convex set X and 
X £ X, Consider the two programs % 

Program I maximize f(x), subject to x £ X. 

Program II maximize vf(x*)x, subject to x G X, 

Then x" is an optimal solution of program I iff 
X is an optimal solution of program II, 

Linearization technique for solving such non-linear 
programmes attempts to solve a sequence of linear programmes 
in order to arrive at the required linear programme, Bhatt [26] 
has given a linearization technique for solving pseudo- concave 
programmes* 
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In this chapter, ws consider pseud o-montonic 
proararames of the type 

max f(x) 

subject to Ax = b 
(P) 

X > 0 
X e Q d 

where A is an mxn matrix of full row rank, b is m-vector 
and f is pseudo- monotonic function on Q, ¥e also assiime 
that '^xj Ax = b, x> 0^ is bounded. The set Q is such that 
when f(x) is linear, the solution procedure for solving the 
problem is available. For example, the set Q may be set of 
integer solutions or zero - one solutions, 

Kortanek and Evans theorem is not valid if programme I 
is programme P. In view of this, we prove a theorem for 
program ? similar to Kortanek and Evans theorem, in section 2,2, 
Using this theorem, a linearization algorithm is given in 
section 2,2 to solve pseudo-monotonic programmes of the above 
type. In sections 2,3 to 2,7, the linearization technique 
has been used to solve the following problems with 
pseudo-mono tonic objective functions : all integer and 
mixed integer programmes, zero - one programmes, transportation, 
assignment and traveling salesman problems, set covering and 
set partitioning problems. 

Some research workers who have earlier worked on related 
problems are : Swarup [l4l,143],Grunspan and Thomas [72], 
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Florian and RoMllard [53 ], Taha [147] jSuresh Chandra [28 ], 
Anzai [3] , Ishiy Iharaki and Mine [85] , Arora and Purl [4], 
Prohlsms solved by then are for fractional objective fimctions 
which are special cases of pseudo-monotonic functions [ 116] • 
Thus the linearization technique given here is more general 
and solves a large class of problems. 

2*2 The L ineari zation Algorithm 

The proposed algorithm to solve the pseudo-monotonic 
program (P) is based on the following theorem. 

T heorem 2 i Let f be a pseudo-concave fmction on S where 
S = '{s: £ iPj Ax = b,x>0 }o Q. Consider the programmes 

max f(x) 

(PI) 

subject to X £ S 

max Vf (x’^)x 

(PII) 

subject to X £ S, 

If x"* is a solution of (PII), then x* is a solution of (PI). 
Pro of s X* is a solution of program II 

«> 7 f(x^)x* > Vf(x*)x, ¥ X £ S 

=¥ (x-x*) V f (x"'*') <0, X £ s 

=i> f(x) < f(x'^), X £ S (by pseudo-concavity) 

==^ X* is a solution of program I. 
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Theorem 3s If f is a pseudo-monotonic function on S and 
is an optimal solution of 

max vf(x^)x 

(L.) 

subject to X £ S 

such that x^’’’^ £ ^x^jX^, . . , jX ^^' j i*®* for 

some 3,1 < j < ij then x^ , , , ,x^ are all optimal 

solutions of the program (P), 

Proof t Without loss of generality, let x^’^l solves (L^)» 

1 < k < i. Then Vf > vf (z^)x for all x in £. 

As x^^ G S therefore Vf (x^'^)x^'^^ > Vf (x^)x^ which implies 
'v'f (x^) (x^^"*"^ - x^)> 0. 

Since f is pseudo-convex therefore it implies 

f(xk+l) > f(xk). 

Hence f(x^) < f(x^) < < f(x^) < fCx^"^^). (1) 

Now x^ = maximizes (l‘j_)j for some j, 1 < j < i* 

If 3 = i, that is x^ = x^ then by Theorem 2, x^ solves 
program (P), Ifi< 3 <i then 

Vf (x^)x^ > vf Cx^)x for all x G S, 

Since x^ £ S, therefore 

Vf(x^) (x^*-x^) > 0} 
which implies, f(x3)> f(x^) 


(by pseudo- convexity of f)* 


( 2 ) 

(3) 



34 


From (1) and (3) we get 

f(s;j) = = t»«. = f(x^) (4) 

Since pseudo- convexity implies quasi- convexity 
therefore 

f(xj) = f(x^) =»=•> Vf (x^) (x^-x^) < o, or 

Vf(x^)x^ _< vf(x^)x^, (5) 

Combining (2) and (5) and recalling that x^* = we get 

VfCx^)x^ = vf (6) 

which implies that x^ maximizes (L^), 

Now by Theorem 2, x^ is also an optimal solution 
of (P) and so are the solutions x^, j.i-1 

because of (4). 

Theorem 4s If we further assume that the set S is such 

that y £■ S implies max Vf (y)x exists and the set of 

xes 

optimal solutions is finite then the converse of Theorem 2 
is also valid. 

Proof s Obvious 
Convergence 

The algorithm by solving a sequence of linear programs 
generates a sequence of point of S till one of the points 
generated earlier is repeated and in view of Theorem 4, finite 
convergence is granted. 
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Algor ithm ! 

Using Theorems 2,3 and 4 we now state the 
algorithm as given helow s 

Step 1 i Let he an initial point £ Pp such that 
VfCx^) 0* Set i = 1, Solve the linear program 

max Vf(3p)x 

(L. ) 

subject to X £ S 


Let be a solution of (Lj^). 

St ep 2 : If = x^ £ -^x^jX^, . . .,x^^, then 

x^jx^*'*'^} . . .,x^ are solutions of program (P) (Theorem 3), 
Otheivise, go to step 1 with i = i+1, 

Eejnn]^ : If the problem is of minimization in place of 
maximization the statements and proofs of Theorems lj2|3 
would be similar except that the role of pseudo~convexity 
and pseudo-concavity gets interchanged and thus the 
algorithm would work, as such. 

In the following sections we show applications 
of linearization technique to solve some pseudo-monotonic 
programming problems. The pseudo-monotonic program 


max (min) f(x) 
subject to X £ S 

is referred as pseudo-monotonic integer program, zero-one 
program, transportation and assignment problem, traveling 
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salesman problem, set covering and set partitioning 
problem depending upon the nature of the constraint set. 

For example the following program is referred as 
pseudo-monotonic integer program (all integer or mixed integer); 

max f(x) 
subject to X £ S 

where S = Jx £ jR^{Ax = b, x > 0, x is integer or mixed integer 

“ (as the ease may be) 

where f is pseudo- monotonic on S and S is bounded. 

The fractional programs max(cx*H=()/(dx+8) subject to 

X £ S such that dx+p > 0^ x £ S are special cases of 

pseudo-monotonic programs [116] , At some places we are 

taking examples with fractional objective functions, 

2 , 3 Pseudo- mcaiotonic Integer P r ogramming 

Methods for finding integer solutions to linear 
fractional fimctional programs which are special cases 
of pseudo-monotonic programs are given by Swarup [ 141] j 
Grunspan and Thomas [72] , Anaai [s] , Ishii, Ibaraki and 
Mine [85] . In general, the methods for solving pseudo-monotonic 
integer programs are not available, In[18^9]an algorithm 
for pseudo-monotonic integer programs and some other 
programming problems is given. 
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In this section we use the raethodology developed 
in section 2*2 to solve the pseudo-monotonic integer 
programs. 

Case 1 ; All I nteger- Ps eu do- n i onoton ic Progr am 


E xample 1 


(P) 


max f (x) = 


x^+ 2 -X 2 ■/ (Xj^+2)^ + x| - 1 

O 

(2C^+2) + 


subject to 


x_ 


+ 2x. 


2x, 


X, 


< 3 

< 5 

< 2 


X is all integer > 0. 


Now for any x 


(oc^+l)^ Vf(x) 


"" -(x^+2-X2o() Cx2+(x^+2)o(3 ”1 

- [ 1 ^^+ (3:^+2 )^c“] 

0 

0 

0 

L ° 


^MStos^flief® fimction of this example is due to Bela 
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where = ■/(x^+2)^ + (xp)^ -• 1 

and x^jXgjXg are slack variables. Let S be the solution set. 
Itera tion Is 

Step ! s Let = ( 33030 , 0 , 5 , 2 ) be the initial point, 
bolve the integer linear program 

max Vf(x^)x: = -25/^ Xt-GOOx^ 

CL. ) -L 2 

subject to X G S 

by Gomory’s all integer method C66,743* 

y? - (0,0,1, 1,5,1) is optimal solution of (hj_)« 

St ep 2 i I * 

Iteration 2 s 

Step 1 : Solve the integer linear program s 

max Vf(x^)x = -4/3 Xt - 12x„ 

a„) ^ 2 

subject to X G S 

x^ = (0,0, 0,3, 5, 2) is optimal solution of (L^) 

Step 2 : x^ 

Itera tion. 3 s 

Step 1 j Solve the integer linear program s 

» * * 

3 

max Vf(x )x = -4./^ x-,-12x„ 

(L ) J- ^ 

^ subject to X G S 

which is same as program (Lg). Its optimal solution is 
x^ = (0,0, 0,3, 5, 2). 
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Step 2 t Since = :^ and thus the optimal solution of 
pseudo-monotonic integer program (P) is = (0,0, 0,3, 5, 2) 

-which gives optimal value of f (x) equal to ^ . 

Cas e 2 : Mi xed Integer Pseu do- mon oto njjq. P rogram 

Px^-Xg+l 

max f(x) = 

(P) Xj_+Xg+S 

subject to 

33^+7x2 <20,5 

4X|j^+5x2 < 21 

x^,Xj 3 integers > 0 

For any x, 

(Xg+l) 

- (Xj+1 ) 

0 
0 

where are the slack variables. Take x^=(0,3,0,0) as initial pt. 

Iteration 1 i 

Step ! : Solve the mixed integer program 5 

. max vf(x^)x = 4Xn-Xrs 
(%) ^ 2 

subject to X e S 


(Xi+Xg+2)' 


Vf(x) = 
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by Gomory’s method for mixed integer programs [67j74] 

This gives x^ = (5j0j6,5,l) as optimal solution of (L^) . 

Step 2 : x^ ^ *£^^3^ 

Itera tion 2 s 

S-tSa. i i Solve the mixed integer program i 

? 

max Vf(x")x = x,~6xo 
(L ) ^ ^ 

subject to X e S 

Solution is x^ = (5, 0,5, 5,1) 

St ep 2 ! Since x^ = x^, therefore x:^ is the optimal 

solution of pseudo- mono tonic mixed integer program (P) 
with optimal objective function value = ^ , 

2*4 Ps eudo- monotonic Ze ro -One Programs 

Balas [8] , Geoffrion [60] and Glover [61] 
have given methods for solving linear programs with 
zero-one variables. Florian and Robillard [53] and 
Taha [147] have given methods for solving zero-one 
fractional programs. In general, methods for pseudo-monotonic 
zero-one programs are not known. Linearization technique 
can be used to solve pseudo- mono tonic zero-one programs. 

This solves, as a special case, the linear fractional 
zero-one programs also. 



..Z.ero~.On.e Pr ogram 


Zx^+S^Cg- 4X2” x^+Xg+l 
iillLXl f (X} « ■«» *—» » 

Xj_+2x2+3C2+3x^+2Xg+2 

(P) 

subject to 

x^+ x„ + x„ + 2x. + Xj- <4 
1 2 3 4 5- 

7x^ “■ '^^4 ^ 

11Xt-6x_ + 3x' -3x- 

■i- 2 4 5 

X = o or Ij j = Ij2j*,..j5. 

u 


^ iC* 

T:r Co = lj2, .,,j5), after omission of common positive 
0 

p 

multiplier 1/ (denominator) are as follows % 


ff;L= 

= “X2_+ll3!:2+llx^+4Xg+4 

s 

11^ = -6x1-11x2-11x4-9x5-9 
Ifj = -7x1-11x2+11x3-5x5-6 

= - 3 X 1 - 4 X 2 + 9 X 3 + 6 X 4 

5 

vf(x) = ,..., 4 ^ ) 

dxj^ dx^ 


I V 



Ite ration 1 : 

S tep 1 s = (l,l,0,0j0) is initial point 

Solve the linear zero- one program : 
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min Vf(x^)x = x^“X2-17xg“18x^-7Xg 

(Ll) 

subject to X £ S 

by Balas* additive algorithm for zero-one programs [8,60|148]. 
This gives x^ = (0,0, 1,1,0) as the optimal solution. 

Step 2 ; ^ {x^} 

Ite ration 2 : 

Step 1 : Solve the linear zero-one program t 

min Vf (x^)x = 13 X 2 ^+ 22 X 2 - llx^+llx^+Mx^ 

(L 2 ) 

subject to X £ S. 

We get x^= (0,0, 1,1,0) 

Step 2 i Since sP - x^, therefore x^ is the optimal 

solution of pseudo- monotonic program (P) with optimal 

objective function value = - £ . 

3 

Remark 1 : The initial point x^ = (1,1,0, 0,0) was the worst 
solution with objective function value = £ , 

Exam p le 4 ; 

X.j+2-X^ (Xy^) + X^-l 

Uj+2)^ + 3 ^ 


min f (x) = 
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subject to 


+ Xg < 2 

= 0 or 1 

is a pseudo- mono tonic zero-one program. It can be seen 
that this problem has two optimal solutions, namely (1,1) 
and (0,1) with minimum value of objective function 
being equal to zero* 


2*5 Pse udo- m ono tonic T r ansportation and Assignment Problem s 

Swarup [ 143 ] and Suresh Chandra [28] have given 
methods for solving transportation problems and assignment 
problems with linear fractional objective function. Linearization 
technique can be used to solve transportation and assignment 
problems with pseudo-monotonic objective functions which, 
as a special case, solves fractional transportation and 
assignment problems also. For the purpose of illustration 
we present fractional programs* 


Exampl e 5 s Fract ional Tra nsportation Problem 


(P) 


min f (x) 


3 4 

S 2 c. .X. . 
1=1 j=l ^3 


3 4 ^ 

2 2d. .X. . 

i=l 3=1 
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subject to usual linear transportation constraints with 

availabilities 10,8,6 and demands 7, 8, 5,4, The cost 

matrices /c. and /d. are given below s 
I. ^3/ L ij J 



We take x?- ^ initial basic feasible solution 

X J 

given by north-west corner rule [76 ] , as the initial 
point for the algorithm and is given below : 


(Tableau' 1) 


24 

• 7 

8 

5 

4 

10 

@ 




8 



(D 


6 


1 

® 

(3) 


The circled entries denote the values of basic cells, the 

/ 

first row and first column denote demands and availabilities. 
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Iterat ion 1 s 

b tep 1 ; Solve the linear transportation problem ; 

rain Vf(x^)x 

CL. ) 

subject to X £ S 

where S is the corresponding transportation polytope. 
The cost matrix /vf (x^)^ for the problem is given in 
tableau 2. 


(Tableau 2) 



48 

86 

96 

—66 

^vfCx^)]- = 

10 

-8 

-66 

-18 


-18 

—66 

I 

. 10 

-94 


Dantzig*s uv method [76] gives the optimal 
transportation solution as given in tableau 3, 

(Tableau 3) 


' 24 “ 

■ T‘" 

8 

5 

■“* 4 " 

10 




- 0 - ■ 

8 





6 

1 


© 
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Step 2 : 

\ « mm ni t 

It£ra.ti_m 2 
Step 1 : 


» 

« 

Solve the linear transportation problem : 
min vf(x^)x 


(Lg) 

sub;i6Ct to x G S 


The cost matrix 


{ 


Vf (x^) 


is given in tableau 4. 


(Tableau 4) 



54 

82 

108 

-30 


26 

50 

-30 

24 


: 24 

-30 

26 

-32 


The optimal solution x^ of (L^) is given in tableau 6. 

(Tableau 5) 


24 

7 

8 

5 

4 

10 

© 




8 





6 
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St ep 2 t Since = x^j therefore x^ is the optimal 
solution of pseudo-monotonic transportation problem (P) 

■with optimal value = . 

14 


Exa,nm le__6 


The Fr actiona l As s ignment P ro blem 


(P) 


min f(x) = 


4 4 

5 2 c. .X. . 

i^l _0^]^ ^ ti 

4 4 

2 2 d. .X. . 

i=l j=l 


subject to usual linear assignment constraints, - 

Where the cost matrices given below 
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The extreme points of the assignment polytope of 
order n are the permutation matrices of order n. Corresponding 
to a permutation ir £ the symmetric group of order n, we 
can associate a permutation matrix ZCtt) = hi as follows ; 

X. .(it) = 1 iff irCi) = 3 
10 

= 0 otherwise. 

Furthermore, we mean by the symbol ij^ig...,!^ as the 
1 2 • t • n 

permutation ( . . . ), Henceforth we represent the 

^1 ^2***^n 

assignment solution by a permutation in one row only (base row 
understood in its natural order). 

Iteratio n 1 t 

Step 1 ; Let TT^ = 1234 be the initial assignment. 

Solve the linear assignment problem s 

min Vf (Z(Tr^ ) ) x . . 

(L ) 

suboect to X G S 

where S is the assignment polytope of order 4 whose cost 
matrix is as given in tableau 6, 
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(Tableau 6) 



The Hungarian assignment technique [98] gives = 4123 
as the optimal assignment to (Lj)« 

Ste^2 i TT^ t 
Iteration 2 i 

Step, 1 s Solve the linear assignment problem : 

min VfCXCTT^)) x.. 

-i- J 

(Lg) 

subject to X G S 

■whose cost matrix is given in tableau 7. 

(Tableau 7) 



3 

Optimal assignment to (Lg) is ir = 4123 

step 2 ; Since = ir'^j tlierefore optimal assignment 

solution of pseudo-monotonic assignment problem (P) is 

q c: 

TT = 4123 with optimal objective ficiction value = - . 

6 

Bejiai^^ : The Initial assignment = 1234 with which 

3 

we started was the worst assignment with value 2^ , 


2*6 P s eud orjiionj^^^^ jE r ax® ling^ 


The traveling salesman problem of order n with, 
linear objective function is vrritten as 

n n 

minimize fCx) = 2 S c. .x. . 

i=l j=l 13 13 



subject to 

n 

2 X. . 

j=l 

= 1 

(i = 1,2, .. . ,n) 



n 

2 X, , 

i=l 13 

= 1 

(j — l,2,.*.,n) 



X. , = 0 

13 

or 

1 


such that the solution set of those x.-.’s for which = 1 

13 ij 

forms a tour where a tour is a sequence of integers taken 

from (l,2,,,.,n) in which each of the integers appears 

once and only once and we mean by a tour ir = (i;j_,i 2 , . . . , ip^ . 

the solution X;; . = ^ = ... = i Xi i= 1 

il-ig 12-^3 ipip+1 In 1 

and vice versa. Therefore a tour on n symbols can be 


represented by a cyclic permutation on n symbols and thus 
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tilers are (n-l)l tours which are all feasible solutions of 
the problem (Pq). 

The usual terminology is that the n integers correspond 
to n cities or nodes of a graphj is the distance 

matrix*' giving "distances” from node i to j or the arc 
lengths i to a tour it being a closed path passing through 
each node exactly once and the objective function gives 
the sum of the arc lengths over the arcs included in the tour* 

Traveling salesman problem arises in number of contexts* 
To name a few of them are routisag problems, job scheduling, 
computer wiring[22,l0^ and minimizing wall paper wastage [ 56 ] . 

The proposed traveling salesman problem with 
pseudo-monotonic objective function is 

min f (x) 

(P) 

subject to X G TP^ 

where TP^ is traveling salesman polytope of order n 
which is the convex hull of the cyclic permutation matrices 
and is the convex feasible set of the problem without 
integer restrictions. 

There are several methods, exact as well as 
heuristics, available for solving the traveling salesman 
problem with linear objective function, ¥e essentially 
need the exact methods for our purpose. See [22] for brief 
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discussion of them. In general, the methods for 
pseudo-monotonic traveling salesman problems are not 
available , \Je here use the linearization 

technique to solve such problems. The branch and bound 
method of little et al [l09] was used to solve linearized 
problems. 



min f (x) 


6 6 

2 2 c X + c( 

1=1 j=i ^0 

6 6 

2 2d. .X . + p 

i=l ^=1 ^0 10 


subject to usual traveling salesman constraints such that 
6 6 

_2 _2 dj^^x^^. + p > Q for all tours £ Pg, is an example 

3-*"l ^ ““1 


of 6x6 pseudo-monotonic traveling salesman problem. «< and p 
in the present case are taken to be zero matrices. The 
matrices gi'ven below 
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Ite 1 : 

5-tejB, 1, i Let = (1,2,3 ,4,5,6) be the initial tour. 
Solve the linear traveling salesman problem : 


ap 


rain Vf(x^)x 
subject to X e Pg 


for which the distance matrix is given in 

tableau 8, by branch and bound method of Little et al [l09]. 
The optimal tour to (L^) is x^ = (1,4, 6, 5,3, 2) , 

Step 2 : x^ f 

Iteratio n 2 : 

Step 1 s Solve the traveling salesman problem : 


min vf(x^)x 

(Lg) 

subject to X e P 



for which the distance matrix 


is given in tableau 9, 
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Step ^ : Since therefore y? is the optimal 

tour to pseudo-monotonic traveling salesman prohlem (P) 

which gives optimal objective function value equal to - , 

5 

Re mark 3 s The above example was also tried by subtour 
elimination method, of Eastman [49] by starting with same 
initial solution x^ = (Ij9j3,4,5,6), But, however, the 
first assignment solutions in both iterations applied to 
tableau 8 and 9 gave the tours x^ = (1,4, 6, 5, 3, 2) and = x^. 

2,7 Pseud o-monoto nic Se t_ Cov eri ng and Set Partitioning Problems 

The set covering problem (SC) with linear objective 
function is : 

min z = cx 

subject to 

Ex > e 

and X, = 0 or 1 ( j = l 3 ...,n) 

J 

where E = H} ■is a mxn matrix of zeroes and one's, e 
is m~column vector of all l^s» The usual interpretation is 
that the columns' of E can be thought of as sets and the rows 
of E as elements such that 

e^. = 1 if ith element £ set j 
= 0 otherwise 

The SC problem is to find a cheapest union of sets from E 
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that covers evovy element of e. 

If the inequality constraints are replaced hy equality 
constraints the corresponding problem is set partitioning 
problem (SP), Thus in SP problem, we seek a cheapest union 
of disjoint sets from E which covers e« SC and SP problems 
have many practical applications. For details refer to [58, 135] . 

Balas and Padberg [9] , Bellmore and Ratliff [23] , 

Etcheberry [52] , Garfinkel and Nemhauser [57] , Jenson [ 06 ] , 
Lawler [lOl] , Lemke, Salkin and Spielberg [l07 ] , 

Salkin and Koncal [136,137] have given methods for SC 
and SP with linear objective functions. 

There are some models where the SC and SP problems 
with fractional functions arise. For example .the locational 
models where the objective may be to minimize the ratio of 
cost incurred to some measure of public facility. For 
example the construction of rationing depots in a big city 
where a particular depot is assigned to demands of users 
of particular sectors only. The public utility of the 
depot depends upon the number of persons who find convenient 
to visit it. For details of locational models refer to [130,149] . 

Recently, Arora and Puri [4 ] have given an enumerative 
technique for SC for those fractional objective functions in 
which the cost coefficients c.>s of the numerator are non-negative. 
In general, the methods for pseudo-monotonic SC and SP have 
not been given. The linearization technique is used here to 
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solve SC and SP with pseudo- mono tonic objective functions. 

This solves, as a special case, the SC and SP with linear 
fractional functions also. Furthermore, unlike Arora and 
Puri [ 4 ] , there is no restriction on the cost coefficients 
c^*s, ¥e give below the solutions to two examples, one ea.ch 
for SC and SP. 

Example 8 : Fra ction a l Set Covering Problem 

2x^+X2+3x3+2x^+X5+Xg+PXY+4x0+2Xg+2x^Q 

2x^+2X2+2Xg+3x^+3Xg+3Xg+4x2+4Xg+4Xg+3x2^Q+l 


subject to 


where S = -f x £ R^lEx> e, 
matrix E is given below : 


X is 


zero*»one vector 


> 


and 


1110 10 
0 0 0 0 1 1 
0 0 1 0 0 0 


10 0 0 
0 0 10 
1110 


E = 


10 0 
0 0 1 


0 0 0 


0 0 1 


0 0 
1 0 
1 0 
1 0 


0 0 
0 0 
0 0 
0 0 


0 0 1 
10 1 
10 0 
0 0 1 
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df 


dXj ^ 3 - lj2j,,,,10), after omission of common positive 




multiplier 


w I m i i w i . Ill 

(Dominator 


are as follows s 


af 

a”x^ 

of 

ax„ 


of 

It., 

ox. 


2Xg- 2x3+2x^+4Xg+4Xg+4x^+4Xg+2x^Q+2 


■EXi-4X3-X4+X5+Xg-4Xg-.^j_0+l 


= 2x^+42C2+5x^+7xg+7xg+82;^+4Xg+8Xg+5x^Q+3 
= -2x^+x2-5x2+3Xg+3Xg+2x^-4Xg+2Xg+2 


PX- 

ax5 

SX- 

dXa 


- -4x2_-X2-7x3- 3x^- 23c^-- 8Xg~ 2Xg- Sx^^+l 

_ ^ 2^ , IX. = 2^ , IX = 2IX- , = Xt 

/^TT- * Pv_ /?Y ^ r'-TT J pTr_ * ?^y P^y 


CXy OXg C“Xg 


0x9 -ox^ ax^9 


Iter_ation ' 1 : 

S tep 1 : Let x^ = (1,1, 1,1, 1,1,1, 1,1,1) be initial point 

of (P). Solve the linear SC problem : 


(ip 


min Vf(x^)x = 22xj^~9x2+53s^+2x^-29Xg“29xQ-18x^+44Xg“18x9+2x. 
subject to X 6 S 


10 


where S is the feasible region of SC, using Lemke,Salkin 
and Spielberg^s [ 107,135 ] enumeration procedure. 

The solution is x^ = (0,1, 0,1, 1,1, 1,0, 1,1), 

SteP.X *• ^ “C^ll 
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It eration 2 s 

• Solve 


Vf (x®)2: 





subject to X e S. 


The solution Is = (0,0, 0,1, 1,1, 1,0, 0,1) 

S tep 2 t f 

I terati on 3: 

: Solve tain 'v^(x2)x = 18x^+r^+35x^+10x^>‘7:iz^-7xQ+2^ 

+36xg+2xg+l0x^^ 

(Lo) 

subject to X 6 S, 


The solution is x'^ = (0,0, 0,1, 1,1, 1,0, 0,1) 

4 3 4 

S^' j j ^ep 2 : Since x = x , therefore x is the optimal 

solution of pseudo-monotonic SC (P) with optimal objective 

function value equal to .JL-. 

17 

Rem ark 4 ; In SC problem it is assumed without loss of 
generality that the costs are positive, because non-positive 
costs can be handled by setting x. = 1 for c. < 0 and deleting 
row i such that = 1, Thus, because of the appearance of 
the negative costs in the objective function of the linearized 
problems, practically it is seen that size of the reduced 
problems with positive costs is considerably less and many 
of them get solved just by inspection* 
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Frac tiona l _S_et Partitioning Problem 


min fCx) 
(P) 


X1+2X2+X3+3X 


2x^+X2+2x2+x^+3Xg+4x6+Xp+2Xg+3Xg+x^0'^^ 
subject to X £ S. 
where S = ^x £ R^jEx = e, x is zero-one vector^ and the 


matrix E is 

C" 

i 


E = 




0 0 1 
0 0 0 


0 

1 


1 

1 


0 0 
1 0 


0 10 10 0 0 
0 0 1 1 0 0 0 


0 

0 

1 


0 0 
0 0 


1 0 


11000010 


110 


0 


0 


1 

0 


0 

1 

0 

0 


pX- (j = 1,,..,10), after omission of common positive 


dXi 


multiplier 1/ (Denominator are as follows 


ox. 


• 3 X 2 “ 5x^+Xg+2Xg- 3x^- Xg- 7 x^Q +1 


PX_ = 3Xj+3x2~x^+5Xg+7Xg+3Xg+4Xg-2Xj|^g+2 


5x, 


gf 

gx^ 


gf 


gf 

gX3_ 


^ a* 5xi+X2+5x3+8X5+llXg+X7+5Xg+7x9-Xjo+3 
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C’X 


5’Xg -Sxj-TXg- 2x2-llx^-x:g-7xY~2xg“5xg~15xjQ+l 
x^ o'Xg ^ cxg ox. 


df_ - ^ 
0 


f. ^^ai£*«sa» “*" ^.a^Lr: 

2 “^1 


Bf 

= 7xj_+2x2+7x2+x^+lls:g+15Xg+2s7+7Xg+10Xg+4 


oxio 

I teration 1 : 

Step 1 : Let x^ = ( 0 , 0 , 0 , 0 , 1 , 0 , 1 , 1 , 0 , 0 ) Le initial point of 

(P), Solve the linear SP 


(L^) 


min vf (x^)x = - 2 x 2 _+ 8 x 2 - 2 Xg+ 14 x^- 6 Xg“ 10 Xg+ 8 x;^- 2 Xg+ 20 x^Q 
subject to z £ S. 

I 

by enumeration procedure of Garfinkel and Nemhauser [ 57 ] • 

The solution is x^ = ( 0 , 0 , 1 , 0 , 0 , 1 , 0 , 0 , 1 , 0 ) 

Step 2 s x^ »w 
Iteration 2 : 

St ep 1 : Solve min vf(x^)x = 2 x 2 _+ 16 x^+ 2 x^+ 26 x^» 2 x^'- 6 xQ+ 16 xy 


(Lg) 


+2Xg+8Xg+36x^Q 


subject to X G S, 


The solution is = ( 0 , 0 , 1 , 0 , 0 , 1 , 0 , 0 , 1 , 0 ) , 
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Step 2 : Since therefore is the optimal solution 

of pseudo- monotonic SP problem (P) with optimal objective 
function value equal to “ , 
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SMSlTiyiTY ANALYSIS OF LINEAR COMPLBIENTARITY 
PR OBLB-I HITH A PPLICA TIONS * 

3*1 In t ro duct ion 

The linear Compl^nientarity Problem (LGP) is to 
find a solution of the system of equations 

w = PIz + q 

¥ > 0, z > 0 (1) 

«»«■• *»«• 

T 

w^z = 0 

where w = (w^jWgj • •• j ^ ~ • • *>23^), M is a 

square matrix order n and q is an n- vector. LGP is closely 
connected with’ mathematical programming, game theory, economic 
equilibrium theory and fixed point theory. 

LGP was first studied by Dantzig and Gottle [42 ] 
in 1967 in the field of mathematical programming. They 
proved that necessary optimality conditions for a quadratic 
programming problems given by Kuhn-Tucker [99 ] become 
the LGP and for convex prograiiming problems the Kuhn-Tucker 
conditions are sufficient as well. Also see [37,38 ], 

They also provided ai algorithm for solving LGP which was 
referred to as the -pri ncipal pivoting me thod. 

*Presented in Joint National Meetog of ORSA/TIMS held at 
Los Angles (U.S.A) in Nov. 1978 [20 J , 
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Another algorithm for LCP is Lemke^s algorithm 

which was originally developed hy Lemke and Howson [ 106 3 
for solving LCP corresponding to the problem of computing 
an equilibrium pair of strategies in bimatrix game problems 
which was later extended by Lemke[l04 }to solve general 
linear complementarity problems. See [119,124 1 for 
details of Lemke* s Method. The third method is due to 
Graves [S^l . 

All the above three methods [37,69,104] are 
pivoting algorithms and are also referred to as compl^^menjsary; 
pi voting a lgorithms. Under certain conditions on the matrix 
M, these algorithms solve the problem efficiently. Of the 
three algorithms, Lemke *s algorithm is applicable to the 
broadest class of problems. Cottle and Dantzig [37] in . 

1968, Lemke [lOS] in 1970, Saves [50] in 1971 and Saigal [134] 
in 1972 etc. have considered some conditions on matrices M 
under which Lemke *s algorithm either leads to a solution 
or shows that the problem has no solution. Several types 
of matrices M have been defined for which it can be shown 
that LCP has a solution. For details refer to [l05 ] , 

Lemke* s algorithm has several applications in linear programmir^ 
(Ravindran [ 127 ] ), quadratic programming (Eaves- [51 ] ), 

non co-operative games (Wilson [l54] , Rosenmuller [133]), 

Since Lemke* s algorithm, a lot of combinatorial methods have 
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"been developed by using his idea* For example see Scarf* s [l38 ] 
raethod and its applications [l39,140]. 

The second class of algorithms for the LCP is 
class of cutting- plane type or branch and bound type 
algorithms# ICirchgassner [95] has developed a cutting 
plane algorithm while Jeroslqw [ 87 ] has characterized 
all the valid cutting planes for this class of problems. 

Ibaraki [83,84 ] has applied branch and bound concepbsto 
solve LCP, These algorithms can solve LCP without any 
conditions, but their computational efficiency may be 
questionable. So far more interest has been centred on 
pivoting algorithms and their degree of efficiency is found 
comparable to that of the simplex method for solving an LP 
of a comparable size. The computational experience of 
Ravindran [128] reveals the superiority of Lemke’ s compKtmentary 
pivot method (Lemke^s method in short) over the simplex method. 
Recent studies by Polito [ 125] and Lee [102] have shown the 
superiority of Lemke’ s method for solving convex quadratic 
programming problems as well. Thus LCP can serve as a unified 
theory for studying linear and quadratic programs and bimatrix 
games. Besides these applications, LCP is used to study minimum 
distance problem, problems in the plastic analysis of structures 
[no], in the elastic flexural behaviour of reinforced concrete 
beams, in the free boundary problems for journal bearings, in tiie 
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study of finance models and in several other areas* See 
[37,119,124] for these applications. 

The post optimality analysis and sensitivity 
analysis are fully developed aspects of linear prograiiiming , 

Need for such a analysis often arises in operations research 
problems because the model parameters used thereof are many 
times rough estimates subject to some ranges of variations 
or sometimes as a result of some policy decisions they need 
to be changed and consequences of such changes are to be viewed. 
In this chapter, we discuss the sensitivity analysis of LGP 
through Lemke's complimentary pivot method. The follovring 
changes have been discussed, 

(i) Changes in right hand vector q. 

(ii) Changes in one or more column vectors of M, 

(iii) Simultaneous changes in q and in one or more 
column vectors of M, 

It is observed that changes in cost coefficients and requirement 
vector of linear programs and convex quadratic programs can 
be more conveniently discussed through LCP, 

As an application of this procedure, we have further 
shown that linear fractional functional programs and quadratic 
fractional functional programs can alternatively be solved by 
considering the sensitivity analysis of their corresponding 
LCP^s by considering changes in q only. 
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In section 3,2 we have discussed the sensitivity 
analysis of LCP, Section 3,3 illustrates the procedure by 
an example considering the sensitivity analysis of convex 
q.uadratic programs. Section 3.4 gives the applications of 
this analysis to solve linear fractional functional programs 
and quadratic fractional functional programs. 


3.2 S.ens itivit v An alysis Of Linear 


Sensitivity analysis of LCP is considered by 
considering the changes of the following types, 

(i) Changes in q vector. 

(ii) Changes in one or more column vectors of M, 

(iii) Simultaneous changes in q vector and in one or 
more column vectors of M, 

We discuss each of these changes one by one. 


(i) Changes in q, vector 

The initial tableau form for the LCP(l) is given 

as below ; 


(Tableau 1) 


Basis 

* * * 

Wg. ,, 

^n ^l'* 

« Z 0 

s 


^0 

q 

Aq 

q+ Aq 


1 . .. 

0 . .. 

0 

Xs 

^.,.m3_n 

-1 

% 



• 

• 

* 

m * 

m 

• 

« 

• 

# 

m 

-• 

• 

c 

m « 

* 

• 

• 

m 

• 

• 

w 

0 

« 

1 

• 

0 -m T 
si 

• • 

ss 

“«isn 

• 

-1 

' 

• 

Is 

• 

Aqg 

• 

‘Is'" 

• 

* 

0 

» 

0 

« « 

^ -°nl 

« 

-^ns 

* 

"■'%n 

•• 

-1 



• 



Let the vector q he changed to q=q+ ^q \'rhere 


Aq = 


Aq, 


AQ 



In the initial tableau, the vector of changes AQ can be vie^#®d 

0-S 



1. . .0. . .0 

• * o ; 


« 

Aq = 

♦ * • 

0 1 0 

1 • • • 


Aq 


: « « « 

0 0 1 


<1 


i.e, ¥e view as if the changes aQ^?.*.} Aq„j.».} Aq„ 

J, w 11 

have been introduced as variables in the initial tableau 
vrith their corresponding activity vectors as e^, . . . ,6^, . , ,,ej^ 
though explicitly they need not be introduced in the tableau 
as the identity matrix itself is present in the initial tableu. 
If the same sequence of arithmetic operations of Lemke* s 
pivot algorithm is carried throughout on the vector of changes 
aq then it can be observed from the final tableau 2 that Aq 

5jC 

traiisforttis to /\q as below ; 
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Aq 


r 




>A. 

‘T’ 


r 

l> 


'^11* ' 

•Pis*. 

* Pin 



« 



* 

s 


« 

• 


« 

* 

« 


* 

sjc: 


P'sl 

^ss 

p’’" 

s-’sn 



• 

• 


• 

«) 

« 

« 

« 

o 


* 

• 

<1 


?;i 

3* 

•°ns 

Pi?n 


1> 


X 




!. 

J 


n 


or Aq„ = 


2 . Aq. , s = 1, . ..,n. 

j=l J 


( 2 ) 


■where is a basis inverse present in the final LCP 


tableau 2, 


(Tableau 2) 


[sasis 

Wl*' 

» • • • 

..Wn 

•^1** 


.Zn 

^o 

>lc 

q 

* 

Aq 

— ^ 
q= q+Aq 

B-^l 

O 

^11 

®!s 

Pin 

“^11 

* 

>1^ 

’^^In 


'p 

ii 

sK 

Aq^^ 

5}C . >1^ 

qi+ Aqi 

• 







non- 

Basic 




B.Vg 

• 

3* 

* si 

St: 

8^c. 

■ ss 

e 

sn 

ill - 

sl 


“^sn 


* 

qs 

>!< 

Aq. 

kA 

q^ Aq* 


Pnl 

pSs 

Pnn 

>lc 

%Ll 

={< 

%s 

=4: 

'^nn 


❖ 

qn 

f 

Aq.* 

qn+ ^ In 


and q = q Aq which is same as obtained by performing the 
same sequence of pivot operations on q. (A q* is easy to 
calculate than directly calculating q because many components 
of Aq may be zero.) 
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Thus, after solving LCP, when we find, that there 
are some changes in elements of vector q, denoted hy Aq, 
we calculate Aq* hy (2) and then q = q’*' + Aq** Now there 
are two possihilities ; 

Cas e & t q > 0. 

The LCP feasible basic vector w.r.t. q is also 
a LCP feasible basic vector w.r.t. q+ ai- 
Case b s q 0 (i.e.for at least one i, q^ < 0). 
Continue the pivot operations of Lemkev s raethod by taking 
current tableau as initial tableau for the changed problem 
till we get complimentary feasible solution, 
p ara metriz ation of vec tor q 

Instead of considering discrete changes in q 
let q vary continuously as a linear function of some 
parameter ^ i.e. 


^ = q + "^Cf 


where f is some fixed vector, < is scalar, q being 

•* 1 ^ ^ 

changed vector of q. If B is the basis 

corresponding to = 0, then 


B 


J = B'^CqW ) 


= B*"^q + =<B”^f 


q = q 4- 


* ■ oCB~lf 


or 


q = 2 B .f 

s s S 3 3 


the 

inverse. 
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Select 1J.2_ and p,2 such that 



- oo 


if for all s 


for those s for which 

2 6^.f. > 0 

3=1 3 

2 ^*.f < 0 

3=1 3 ~ 




min 

s 


oo 


-1 

n 

2 S . 

3=1 3 

if for all s 


t 

7 Y 

fn-^f-- J 


for those s for which 

2 < 0 
3=1 ^ 

2 P* f . > 0 • 

3=1 ^ ^ 


For < oc < p.2,th6 final LCP feasible basis w.r.t. q is 
also a feasible basis w.r.t, q4c<f which is the case a above. 
For ®C [ p,2_} j we have case b . The sensitivity analysis 
with regard to both cases has been studied above, 

(ii) Changes in colun'ins of M 


For the sake of simplicity we first consider changes in 
one column of M. Suppose that the kth column of M is changed. 

Let the change be 


A 




* 

^“sk 
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This change, in the initial tableau cau be viewed as follows i 



i,e. it may be viewed as if are 

vari8.bles with their corresponding activity vectors as 
®1» " '®s’ • ' ’®n? 'though we really don*t have to introduce 
them in the tableau as the identity matrix; itself is present 
there. 

By carrying out the same sequence of pivot 
operations of Lemke’s method on the vector of changes Am^^ 
it can be observed from the final LCP tableau 2 that 



or 


= 2 Bot 

3=1 ^ 




and 


“k + ^“k 





s = 1,2, ... ,n (3) 


After affecting this change in final LCP tableau, we 


have two possibilities 
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fiase .a ; If kth column is such that the corresponding 
variable is a non-basic variable then the final complimentary 
basic feasible solution w.r.t, old matrix continues to be 
complftmentary basic feasible solution w.r.t, changed matrix. 
In fact in such a case there is no need to calculate m|^, 

* Now suppose that the kth column is such that 
the corresponding variable is a basic variable in the 
final LCP tableau. If kth column is in t-th place j say, 
of the basis then before changes 

= 1 

mlk = 0, i ^ t 

After changes, the kth column looks like 

% = 

This effects the basic structure of the tableau. To reduce the 
tableau to its proper form, we pivot the current tableau with 
respect to (t,k)th element where by a proper form we mean 
the form of the tableau where a basic vector has a unit 
element corresponding to its place of appearance in the 
basis and zeroes elsewhere. However, we have assumed above 


Am! 


sj: 

Ik 


1 + A m. 


tk 
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tliat We can always pivot (t,k)t]i element i.e. we have 
assumed that 



- 1 . 


The results of pivot operations on q" are 


^i ” 


A * 

(1+ ai4,j) 





(1+ Am*^) 


where q is the transformed vector. If q > 0, the current 

LCP basis continues to be complementary feasible basis. If 
q |l 0, restart the pivot operations of Lemke^s method 
taldLng final tableau as initial tableau for the changed 
problem till we get final complimentary basic feasible 
vector* 


Changes in more than one column of M can be 
treated the same way. First we affect all the changes in 
the LCP tableau and then transform the tableau to its 
proper form. Finally we deal with the transformed q* 


(iii)_^Simul^ane^ cjl^^e.s_iA and^ 

Changes in columns of M and vector q can be 
considered simultaneously as well. At first, changes with 
respect to all columns of M and column q should be affected 
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together in the final LCP tableau. The pivot operations 
are applied to deal with the disturbed tableau to restore 
it in the proper form corresponding to the positions of 
its basis vectors and then finally treat the q vector, \Ie 
again assume here that pivoting is possible at each step* 

3,3 

Sensitivity Analysis of Quadrati c Programming Prob lem 
Boot [27 ] has considered the sensitivity analysis 
of convex quadratic programs by considering infinitesimal 
small changes in the parameters of the problem. However 
the results are approximately valid for small changes of 
the problem parameters. Me consider here the post-optimization 
of quadratic programming problems by first converting them to 
LCP, The analysis takes cape of all types of changes, small 
or large. The procedure is illustrated below ; 

E xample 1 i Consider the convex quadratic program 

min f (x) = ~6 x^+2x|-2x^X2+2x2 

(I) 

subject to 

-X]_-X2^ -2 

x_ ^ 0, x^ > 0 

1 ^ 2 •• 


Viewing the objective function and the constraints in the 
standard form 
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1 T 

min fix) = cx + ” x Gx 

2 

subject to 

Ax > b 

0 

the corresponding LCP is derived below : 

-6 

0 
2 

and LCP is w-Mz = q 

The tabular foira for which is (tableau 3) 



— 


G -A^ 


4-2 1 . 


c 




-2 4 1 

} G = 



A 0 


-1 -1 0 





i 



( 


(Tableau 3) 


Basis 

Wi 

^^2 


^1 

^2 

^3 

Z 

0 

q 

Aq 

Wi 

1 

0 

0 

-4 

2 

-1 

-1 

-6 

^4 

■ 

Wg 

0 

1 

0 

2 

-4 

-1 

-1 

0 

F 

i 6 

, Wq 

0 

0 

0 

1 

1 

1 

0 

-1 

2 

4 


( ^q column shom here for convenient future calculations) 

Solving by Lemke*s Complementary Pivot Method [ 119 ], the 
final LCP tableau is (tableau 4) 
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(T clToI© 3A1 4 ) 



Change in a ; Suppose as a result of some policy decisions 
the quadratic program is changed as 

min f (x) = -10x^+6x 

(II) subject to 


g+2xj- 2x^Xg+2x| 


'^r^2 


l’ 2 


> -6 
> 0 


Here M remains the same and 



"-10 " 



" -4 

q + Aq = 

6 

so that 

Aq = 

6 


6 



4 

U. 


To avoid space and for the sake of computational convenience, ai 
column is listed in tableau 3 as sho-wn. This is particularily 
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convenient if the LGP solutions are needed for set of values 
of because in that case pivot operations may 

simultaneously be carried over these changes. IJe have 



Regarding this final LGP tableau (4) as initial tableau 
and applying Lemke*s procedure to affected tableau we 
reach the following tableau 5, 


(Tableau 5) 
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Optimal solution to changed primal problem (II) is 


* - 5 . = 0. f(x^') = -25 


xi = 


2*2 . ■ 2 

It took three pivot opera.tions by post- optimization 


procedure to solve the changed problem and two steps when 
the problem was solved all over again witli the changed datai 
However 3 in this casej the change a 1 '^■fas a very large 
change with respect to q. 

Let us now consider changes in elements of M and 
q both after we have solved program (I), 

Consider the changed quadratic program 


min f(x) = -6 xj+2x^-2XjX2+2x| 


(III) 


subject to 


-x^-2x2> -4 


X, ,x^ > 0 

JL 2 - 


Therefore 



4 

-2 

1 


-6 

AM == - 

-2 

4 

2 

, q+ Aq = 

0 


-1 

-2 

0 

i 

4 



r -t* 

0 


0 


” o~ 

1 


'o’ 

ABj^ = 

0 


0 

, Am^j^ 

-1 

fi 

0 

1 

0 


1 


0 

L J 


2 

u 
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Again for the sake of computational convenience, 
the changes and Am^ may he introduced in the 

initial tableau. We, however, give below the final 
LCP tableau (6) (with regard to matrix M and vector q) 
with changes introducer as below, 

(Tableau 6) 



After affecting these changes we get tableau 7' 

(Tableau 7) 


Basis 







- 

^1 

Wg 

""3 

^1 

^2 

^3 

q = q + Aq 

^3 

1 

2 

2 

~1 

0 

-1 

3 

2 

-1 

^1 

- 

1 

X 

1 



5 

“ 12 

12 

2 


2 

12 

2 

Zp 


. 3^ 

1 

0 

3 

1 

3 

1 , 

1 12 

12 

2 


2 

12 

2 
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Basic vectors in colurims corresponding to Zg and are 
not appearing in their proper form, lo Bring them in the 
proper form, we pivot around the elements (SjS) and (1}6) 
and reach tableau 8. 


(Tableau 8) 


Basis 

*^2^ ^3 ^3 

q 

^3 

2 0 0 1 
^ 14 7 

0 

^1 

> 1 1 „ 2 10 0 

7 14 7 

2 


^ 0 10 

1 

Cj 

14 28 14 



Since q ,> 0, the final LCP tableau is reached. 

tmm 

Optimal solution to primal problem (III) is 
= 2, scg = 1 with optimal value = -6. 

The post~optimization procedure gave the solution 
to changed quadratic program (III) in just two pivot 
operations, needed to restore the tableau to its proper 
form. When the changed program (III) was solved all over 
again, it took three pivot operations (Sometimes it is 
useful in sensitivity analysis to know the set of values 
for model parameters for which the solution remains optimal). 
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Con venient Changes in Quadratic. Programs 

The LCP formulation of a quadratic program indicates 
that changes in hjC and a column of G can he conveniently 
carried out by the method described above. However, 
changes in A are more cumbersome, in that if one column 
of A is changed, then many coluimis of the M matrix are 
affected. Similarly, the LCP formulation of a LP indicates 
that changes in the cost vector c and requirement vector b 
can be conveniently carried out while the changes in A are 
cumbersome • 

3.4 A d plica tions 

3.4.1 Li near Fractiona l Functio nal Progra mming 

Lineavr Fractional Functional Program (LFFP) is 

T 

min fCx) = 

d x+B 

(LFFP) subject to 

Ax > b 

X > 0 

with the condition that d^x+P > 0 for all feasible x« 

Chames and Cooper [si] , Dinkelbach [44] and Martos [ll4l 
have given methods for LFFP. ¥e use the technique of 
Dinkelbach [^] to solve the LFFP. The parameterised 



LPs requix-sd to be solved are handled by considering the 
sensitivity analysis of q in the corresponding LCP. By 
Dinkelhach^s technique we. hc.ve to solve 

min 9 (x, z ) = (c'^x+c( )- z (d^x+8 ) 

= (c-zd)^x+(c(-pz) 


(PIP) subject to Ax > b 

X > 0 

for different values of z till cpixjz) = 0. 
The corresponding LCP is 


w-Mz = q 


where q = 


(c-zd) 

-b 


¥e denote this by LCP(q). 


c 

■b 


Algori thm t The algorithm is described in the following 
steps. 

Step 1 i Set z = z° = 0 and solve LCP(q^) with q^ = 

Let x° be optimal solution to corresponding primal 
problem (PP), Set i = 0. 

Step 2 i (i) If 9 (x^,zi) - q then x^ is the optimal solution 
to (LFFP). Stop. 


(ii) If 


9 (x^,z^) < 6 j where 6 is some pre-as signed 
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positive number then is the approximate solution to 
(LFFP)» Stop, Otherwise, compute, 

^i+l _ _c ^x5 +p{. _ 

d^x^+B 

and go to step 3 with i = i+1, 

S t_ej3_^ _3 i Solve LCP(q^) for 

(c-z“d) 

q = i 

1 I -t) 

by sensitivity analysis applied to LCP(q^_^^) talcing 
vector of changes Aq^^-j = q^-q- t Let x^ be the optimal 
solution to corresponding primal problem, go to step 2* 

Examp le 2 

Below we illustrate the above procedure for (LFFP) 
by an example . 


min f (x ) 


2£,'^2r+T 

1 2 


subject to 

(LFFP) 


Let S be the feasible set 



Iteration 0 


Sis£,l. ! Set z = z° = 0 


(PlPo) 


min 9(x,z°) = 


subject to X £ S 


Corresponding LCP is 

w-Mz = 


where 



r- 



- 




0 

0 

1 

-1 


1 


0 

0 

1 

-1 


1 

M = 

-1 

-1 

0 

0 

»^o = 

2 


1 

1 

0 

0 


-1 


- 





1 


The final LCP tableau (9) for LCP(q^) is 


(Tableau 9) 
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Optimal solution to 

(PLPo 

) is x° = 

(ZijZg 

Step 2 ; 9 (x°,z°) 

^ 0, 

1 1 
z = •= 

4 


Xt e i^aJnjDn 1 




.S.te,p, ,3 : S olve 

min 

9(x,z^) = 

2 




1 

A 


(PLP^) 


(1-7 = qi + Aq 

X O o 


subject to X £ S 


1 

'2‘ 

1 

4 

2 


-1 


Therefore Aq = 


Aq 


1 


4 


- 3 


4 

* Jit =*< 

0 

, q_ = q + Aq^ = 

* ^1 0 o 

0 


1_ J 



1 

- 5 

- 3 

4 

0 

0 


1 

4 

1 

4 

1 


Since q^ ^ 0> LCP tableau (9) continues to be final LCP 
tableau for LCP(q^) also. Optimal solution to (PLPj) is 
= ( 0 , 1 ). 

S tep 2 : 9(x",z^) = 0. Therefore = (0,1) solves 

(LFFP) and its optimal value is = — , 

4 


3,4,2 Quadratic Fractional Functi o nal Programmin g 


Quadratic Fractional Functional Program (QFFP) is 
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1 s:%x+c^z+oc 

2 ^ 

rain — * 

(QFFP) d^x+8 

subject to 


Ax > b 
X > 0 

xiith the conditions that d^x+8 > 0 for all feasible x 

and G is positive semi- definite (PSD) matrix* ¥e use 
the technique of Dinkelbach to solve QFFP, The parametrized 
quadratic programs are handled by considering the sensitivity 
analysis of q in the corresponding LCP, By Dikelbach>s 
method we have to solve 


min (p(x,z) = “ x%x + (c-zd)^x + (cc~^z) 
subject to 

(PCiP) Ax > b 

X > 0 


for different values of z till 9(x, z) = 0, The 
corresponding LCP is w-Mz = q where 


q = 


(c-zd)^ 

-b 


¥e will denote this by LCP(q), 
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Tlie stepa of tlis algorithm to solve this program 

are essentially the same as described in section 3,4*2, 

Exar aple 3 

In the following we illustrate the procedure by an 
example , 


min f(x) 


(QFFP) 


subject to 


2 P 

2x_-2x^x„+2x -6x^ 
1 12 2 1 




*-•2 


0 


Let S be the feasible set of the problem* 
Iteration 0 

Step 1 : Set z = = 0 


min 9 (x, z°) = 


P 2 

2x£- 2 x^X 2 + 2x2“ 6x^ 


(PePo) 


subj< 


to X £ S 


M = 

f ' 

. 

1 

9^ 

4 -2 

-2 4 

1 

1 

f 4n~ 

-6 ’ 

0 


A 0 ; 


‘-1 -1 

0 

0 

2 


Solve LCP(q^) which is the same as that of convex quadratic 
program of section 3.3 and the tableau 4 is the final LCP tableau* 
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Optimal solution to (PaP„) Is =( ^ 1 ). 

^ 2 ^ 2 ' 


s 9(xO,z°) 
I^teratlon 1 


11 


^ 0 , 


= 


U 

6 


(PBP^) 


Solve min 9(XjZ^) = 2x^-2x^x^+2x^+ (- ^9^'^ 


1 

subject to X e S 


'll = 


’o'" 


L25 

} 

1 


11 

1 ^ 

; 


6“ 

11 

, Therefore Aq 


11 

i 6' 

0 


6 

2 



0 

~ 11 ~ 


r 

5 ~ 

6 



6 

0 

« * * , * 

> 1-1 - 1- + m = 
-1- 0 0 


0 

t-J 

2 

0 



1 

_ 



2 


Aq 


Solve LCP(q^) by applying sensitivity analysis procedure 

bo final LCP tableau of LCP(q ), ihe final tableau for 

o 

LCP(q^) is 


(Tableau 10) 


Basis 


^2 

"^3 


"2 

z 

3 

* 

^1 

All 

5|C ^ 

dg = qp Adi 

^3 

1 

2 

1 

2 

1 

0 

0 

-1 

5 

6 

1 

6 

1 

^1 

-1 

3 

-1 

6 

0 

1 

0 

1 

2 

13 

12 

- 1 

12 

1 

^2 

_ 1 

6 

- 1 

3 

0 

0 

1 

1 

2 

i 

1 1 

1 -12 

" fe" 

0 


(Notation ^ after the result of pivot operations leading to 
final LCP tableau, is again referred as ot thouf?h it is its 
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Optimal solution to (PfiP^ ) is M 

J-O 

StejD ^ : 9(x^jZ^) = - .?.§ 0 • Compute c; -2 

72 156 

Iteration 2 

• Solve min 9(XjZ^) = ( 2 x^- 2 x 2 ^x^+ 2 x|}+ (-4Xj+2x2)+2 


(PSPg) subject to X 6 S 



■-4 ■■ 


— 

1 




6 

qg = q^+ Aq^ = 

2 

, Therefore Aqj_ = 

1 

6 


2 


0 


1 


1 

*6 



1 

? 42 " + Aq. = 

1 

" 12 : 

JL 

1 


0 

" 12 




Since > o> tableau 10 is final LCP tableau for 

LCPCq^)* The optimal solution to (PLPg) is x^ = (1,0). 
2 o 2 

step 2 : 9(x ,z"^) = 0. Therefore x = (1,0) solves 

(QFFP) and its optimal value is = -2, 


3.4.3 


Fractional Functional Progr; 


Consider another Quadratic Fra.ctlonal Functional 


Program 
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■g x'^GX+c'^XH-cC 

min f (x ) = — — - - . -*-- - 

A x'^ffic+d'^x+B 
2 


subject to 


Ax > b 
X > 0 

■with 'the conditions that nuraerator > 0, denominator > 0 

for every feasible x, G is PSD and Ji is ND, The problem can 
be solved b^r sensitivity analysis of LCP considering changes 
in q only -which are given by parameterization technique of 
Dinkolbach, We have to solve 

rain 9 (Xjz) = x^ (G~zH)x+(c-zd)^x+Cc<-8z} 

2 

subject to 

Ax > b 

X > 0 

for different values of z till (p(x,z) = 0. 

Co nclusion 

We have demonstrated that sensitivity analysis of 
all those problems which can be formulated as LCP can be 
conveniently discussed through sensitivity analysis of LCP, 
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Furthermore 5 yome of the problems find th'.ir easy solutions 
through changed LCP. A simple computational scheme for 
such problems has been presented. 



CHAPTER - 4 


ML I OJi APPROACH FOR, BONDED .X4RIMM 
LIIiSAR PEOGPJlI''IS* 


4.1 Introd uction 

Tlie Bounded Variable Linear Program (BVLP) is 
max cx 


P(J) 


subject to 


Ax = b 


X > 0 


0 < X. < u., j £ jr 

*“ J J 

where A is an raxn matrix, b is m-vector, x 6 R^, 

J d = •|"l,2, , the index set of variables. Any 

linear program (LP) can also be viewed as BVLP taldng bounds 
on X. to be very large, 

J 

Though the problem can be solved by carryingout 
ra+n equation system ? it bocomes computationally laborious, 
Dantzig [40] has proposed a method to solve the BVLP by 
restricting the basis size to m only. Refer to [76', 119, 148] 
for good expositions of this method. This method has the 


*Th€ work contalnod in this chapter was '^rTs'eiiTed^in Itonua^^ 
Meeting of Bharat Ganit Parishad hold at Lucknow (India) in 
April, 1979 L 90,91 J. 
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following drawbacks ; 

(i) Though the criteria for the selection of the vector 
to enter the basis is the same as in simplex method (i,e, 
the vector a- o.ssociated with most negative z^-c . enters 
the basis)} the criteria for deciding which vector should 
leave the basis is pretty involved. 

(ii) There exist many physica.l situations where the 
optimal LP solution (without bounds) violates only a few 
upper bound constraints. The existing method does not give 
us any information about this, (From here onwards LP 
without bounds is simply referred as LP), 

lie propose tifo methods to solve the BYLP which 
take care of the above drawbacks, to a certain extent* 

The proposed methods essentially fall under the relaxation 
appr oach, i,e, we relax the upper bound constraints and 
then introduce them sequentially. The first method makes 
use of the Dantzig’s method with one upper bound constraint 
onlj^ and the second method is completely independent of it. 

In section 4,2 we prove a theorem which contributes 
to the development of the proposed methods. In section 4*3, 
a method making use of DantzigJs method is given. In section 4,4, 
a dual simplex approach is proposed to solve the relaxed problem. 
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4.2 Main Result 


The BVLP with inde:: set J for the bounded 
variables will be denoted by P(J). We prove a theorem 
below which contributes to the development of the methods 
proposed in the following sections, 

Xheoipeji i If fhe optimal solution to P(cp ) violates 
only one of the upper bound cons traints j say^for 3 = p? 
then in every optimal solution to P(J)) Xp = Up. 

£?j?Pf • Let x° = (xJjXg, bo the optimal solution 

of P((p ) such that x°>Up. Let S(.J) be the constraint set 


of P(J). 

Thus S (<;p>) = lAx = b, X > 0, Xp < Up^. 

For t < Upj define S(p,t) = ^xjAx = b, x > 0, x^ < t 

Clearly S(pjt) O S(p) and S(p,Up) = S(p)* 

Let z,„^^(t) = ma^i: cx. 

xes(pjt) 




Since S(p,t^) C SCpjtg) for t^ < tg < Up, 
it follox^s that z,.,^^(t) is an increasing function of t 

UlcwC 

in the, interval 0 < t < 


5 


> ^ > “p> 


Also because 

z„,„^(t) is a strictly increasing function. Hence x 
in every optimal solution of P(J). 


4,3 _^he_ Jlr tijCl^^^ V_a.r _i pblp Technique 


In this section, we describe an iterative technique 
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to solve P(J) starting withPCcp), LetP(I) denote the 
BVLP 

ina:sc cx 


subject to 

PCD 

Ax = b 

0 < X.. < u. ¥ j e I e j 
j *“ 3 

X > 0 

Let x°(I) denote the optimal solution of P(I) and x^Cl) 
denote the jth component of x°(I), If x9(I) < u. ¥ 3 S Jj 
then x°Cl) solves P(J). Otherwise there exists a p £ J-I 
such that Xp(I) > Up. The essential step of the iterative 
scheme is to solve the B7LP 


max cx 


subject to 

P (I U^P) ) AX = b 

0 < 2 :. < u. •¥ j £ I'|.i{p} 

X > 0 

starting with the solution of PCI), which is achieved through 
the following stops referred to as Substitu tion S teps . 



97 


Substituoion Steps 

cs. JHT rnmtm ii iai ii<i nf!»Tjw 


Substitute X at its upper bound u (Previous 


Theorem). Introduce an artificial variable x 

P 3. 

corresponding to x^ with -H price, M being verj'' large 


positive number at a level = ^p“^p consider it 

a basic variable in place of The column of yji_j’s 

under the artificial variable Xpg_ would be the same as 
that of Xp, 

Recalculate new(z.-c.)’s rovr and the objective 
function value z. 


z = z + (2^-u ) ( -c ) + (2C?-U ) (“M) 

P P y r' P 

= z - Cx°“U ) (H + o) 


p p' 


P' 


If X was present in ith place of basis (the place taken 
P 

by x^„ now) then 

P 


1 


(z.-o.) - Kyij-Vij 


(zj-cp - (M+Cp)y^^, for all j 


^pa~'^pa ^ » 


for the artificial variable. 


Thus Zn~c = - (M+c ) < 0. 

P P P 

Replace x^ by -x' where 0 < x^ < u (The variables x and 31* are 
P P -P-p p p 


Hie t7£Ti" ITonf ihe our"’d~i Sells's ion 'to original 'simp lelT me'&odT 76 X 
only . 
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2cg-Up = (zJ-Uj;) (1) 

* Obtain O ^>'^6 loiow by the 

theorem proved in section 4.2 that x^{l U{p} ) = u . 

Set I = {3 e J|X°(IU (p} ) < M^y. 

• If I = J j z*^(I = x°(.j) otherwise go to 

step 3. 

Con ve rgence 

If ill the optima.1 ‘■;olution to mxn linear program 
without bounds 3 r ba.sic variables are within their bounds j 
then ab che most (in worst case) im-r iterations are needed. 

Thus tho meuhod surely converges. It may become computationally 
laborious in cases where a large number of variables in 
solution X ( T ) exceed their upper bounds. However 3 it is 
conceivable that lH is largOj but iJ-ll,wlth I in step 1, 
is small. In such situations 3 the technique will be useful, 

Iho method is illustra.ted below by an example, 
ma3c z = x^ +4x_+3x^+x +wXr 

subject to 

X +2x +2; +x_^+x^ < 24 

1 2 3 4 5 - 

2Xj+ Xg+dXg+x^ < 48 

x^+3x2+2x2+3x^+Xg< 36 

0 < 8, 0 < x^ < 6, 0 < 2:3 < 8, 0 < x^ < 4, 0 < Xg < 10. 
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_Ite ration ^ 

St ep 1 : The optimal simplex tableau is given below. 


(Tableau 1) 


. -3 


'=3=3 


02-4 



1 


Xn 

J 

x.=12 

5 1 

i 

5 

3 

II 

DO 

2 

3 

^=0 

2 

- 

54 

5 

6 


4 


"1 ^2 


0 

0 

1 

0 


3 


0 

1 

0 

0 


1J5 

3 

1 

3 

7 

3 

7 

3 


T/2^ 

•\r 

o 

.1 

0 

0 


0 


0 

x7 

10 

3 

1 

o 

o 

4 
3 

2 

3 


0 

1 

3 

1 

3 

1 

9 


1 

6 


0 


X, 


8 


7 


3 


i 

3 

4 
3 


5 

6 


Step .. 2 : I = {1,2,4} dJ = ^1,2,3,45^^ 

0:tei)„3 : p = 3, by (1) 

^teji_4. i Substitution steps give the following tableau. 


(Tableau 2) 





'"■T" 

*Q* 

.. , . .. ----- — 




^2 

X3 

1 x^ X5 Fg x^ Fg ■'X3^ 






c =- 3 - 
5 2 

Xc.=12 

5 

1 ° 

0 

. -Ig 1 10 . . 1.7 0 

3 333 

Cq^=~M 

0 ci 



-1 

-i 0 1 1 -1 1 

3 3 3 3 

^2 = ^ 

^3=0 : 

-S 1 

‘ r-v 

1 t 3 

0 

1 0 -1 1 0 ^ 


42 - 4 ^ 


^+ 3 ) 
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(Doulle bar vertical lines are 
from other columns). 

The problem is treated 
only (The notation used is of 

01 = 6 02 = 9 j 

0 = min (6,9j8) = 


put to distinguish the column p 


as BVLP in Vcariables 

Taha[ 148] ). 


Uj^ = 8 


6 


3 


¥e get the tableau 3 (after carrying out usual simplex 
iteration) as 


(Tableau 3) 



©2^ — 12 } 9 2 ~ ^ 3 ^7 — ^ y 0 — 12 


Che next tableau is 


j f-J coll-' cofw 
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(Tableau 4) 


c . 


'IT 

4 i 


’''' 372 '" 

''"”0 "" 


' ‘ o“~'^ 




; 






•^r 


4 

' ■ 1 

•\r 

■“'•4 

' X " 

0 

6 




..5-0 

1 

0 

i 

1 ; 

-3 

1 

3 

0 

-2 


X =12 
/ 

2 

0 

p-3 i 

-1 

0 

1 

1 

-1 


2 

0 

1 

ri i 

l i 

2 

0 

-1 

0 

1 


52 

1 

0 i 

! r 

1 1 ! 

, , . j 

5 

0 

1 

0 

1 



2 

1 

1 

; - 

1 2 


2 



Optimality is readied, 

^ s All the variables are witbiii their bounds. 

So optimal solution is 20 , = 4,:: =0,2c =8jX,=x =0,2=52, 

4 *4 ReXax:ecT,_ Dpial jS3.j;Ty;)lex . Appirjojich. 

In this sectiony we introduce a method to 
solve the BVLP which is completely independent of the Dantzig’s 
method, unlike the one proposed in section 4,3, The method 
relaxes both the lower end upper boimds, as and when required, 
and the dual simplex algorithm is used at all the iterations, 
e^ccept the first iteration where simpl'-x algorithm is used. 

Let x° be an optimal solution of the LP. If 
p G J find x° > u^, then we give below a procedure which 

introduces the constraint Xp < u^. 
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The constraint 0 < z < u is equivalent to 
+ 2 ;^ = Up and 0 5 5 ^p* I^elaxing the upper bound 

on is equivalent to relaxing the lower bound on x^ 
and vice versa. 


If x; = Xp in the final simplex tableau 


of LP then 


X + 2 y. . X. 

P jes J 


where E is the index set of non~basic viiriables. Using 


X +x’ = u„, in •’die above equation, we get 

p p P ? o 


U -Xp = x’+ 2 (-yiPx 

P ^i P jGR 3 


where u -x^ 


u “X^ < 0, 
P P 


Thus X < u can be incorporated in the tableau by affecting 

p ■“ p 

the following changes : 

(i) replacing x° b^^ u -x^ 

P ^ P p 

(ii) replacing variable Xp by x^ 

(iii) changing c to -c 

P P 

(iv) changing sign of all elements in the ith 

row except the sign of ”1'^ in (i,p)th position. 

Ue can observe that sign of (z--c.)*s remains unchanged. 

J J 

The above siteps are referred to as Subs'oitution Stens. 

The resulting problem can be solved by using 
the dual simplex algorithm. The solution obtained will be 
of the problem 
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max cx 

(P) subject to 

Ax = b 
X < u 

p - p 

2: 0 5 a p* 

We observe that Zr, < z. ( by dual simplex iteration) 

where is the optimal value of LP and Zg is the optimal 

value of the program (P)* 

The variable x^ is treated now onwards as a 

bounded variable with upper bound Up in order to impose 

the lower bound on x . 

P 

In general, at any stage of the algorithm, 
the relaxed problem being considered is 

max c:-; 

subject to 

Ax = b 

< u. , j e 

Since this problem is a relaxation of the original BVLP, 
we get 


Zb S 
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where Zg is the optimal value of the original BYLP and 
Zg is the optimal value of the above relaxed problem. 

If the solution to the above relasced problem 
is not feasible for the original BTLPj then one of the 
variables (xj. or violates its upper bound. Such a 
variable is treated by the substitution steps and we 
get a new relaxed problem. 

This new relaxed problem is one of the following : 

k 6 J-Jg 

max cx 

(l) subject to 

Ax = b 

o> j e A-t'iUCk}). 

k G 

max cx 

(II) subject to 

Ax = b 

X. < u. , 

x. > 0- , 

V fum 


} e Ji- (k) 

3 e i„- {k} ) 
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In either case 3 

^ r ’ I' 

where is the optimal value of the new relaxed problem. 
The algorithm can now bo be described as below ? 

AjLgpr^i^m. 

ASPP,.,! i Solve the LPP without bounds by simplex raethod, 

* If* all the variables are within thrir bounds j 

we have solved the BVLP, Otherwise choose p as in (!)• 
■SJig.P.... 3 . : Perform the substitution steps, 

• fipply dual simplex method and return to step 2, 

Since value of the objective function at each 
step is decreasing and the upper bound constraints x- < u.3 

J J 

3 b J which are to be incorporated in the tableau are finite, 
the method converges in finite number of steps# 

finally we point out that the algorithm has 
some drawbacks also, (i) It will take at least one pivot 
operatican of simplex method to remove, the infeasibility 
incroduced by the substitution steps, (ii) If LPP without 
bounds is an unbounded problem then the algorithm may fail 
fco find the so 3 -ution because of its failui'e at the first 
step even though the feasible region of BVLP is bounded. 
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Illustrati 

E xample 


Iteration 0 
Step 1 s 
method [76] 


The method is illustrated by a numerical example, 
max z = 3x^+2Xp 
subject to 

5Xt+2x <70 
1 2 ~ 

“2x^+5x < 30 

1 2 “ 

0 < x^ < 7 

0 < x„ < 7 

- 2 - 

The optimal simplex tableau (5) by simplex 
is given below : 

(Tableau 5) 




3 



0 



X- 

1 


^3 


3 

x^=10 

1 

0 

5 

2 


1 



29 

29 

2 

x„=10 

0 

1 

2 ^ 






29 

29 


50 

0 

0 

19 

4 


29 

29 . 
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.S t ej3_ 2 t Both the variables violate their upper 
bounds and the choice by (1) is bracketed. We arbitrarily 
choose Xj. 


S tep 3 s ¥e (i) replace vf.lus of by = -3 (ii) 

replace variable x-|^ by x^ (iii) change = 3 to c^ = -3 
(iv) change sign of all’ entries in the first row except 
the sign of ”1” in (Ijl) position. The tableau 5 after 
affecting these transformations reduces to tableau 6, 


(Tableau 6) 


c 

% ^ 

/ . 

'-3 

■ ¥ 

o'” 

o"' 








-^B 


'xg' ■ 


^4 

‘1 






~3 

Xt=-3- 

1 

0 

5 

2 ^ 


jl 



29 

*29 

2 

X =10 

0 

1 

2 






29 

29 


50 

0 

0 

19 

4 




29 

29 


Ste p„ 4- : The optimal solution by dual simplex method [76] 

is given in tableau 7, 




109 


(Tableau 7) 


c . ^ 

0 

' ■ Is' ■ 


3^ 

""1 

0 

X3='^ 

29^ 

5 

2 

-44 

2 

5 


193 
‘ s*' 

19 

5 



^ep 2 : The variable x§ = exceeds its upper bound 7, 

d, ^ 

S tep 3 ; Perforraing substitution steps we get the 

tableau 8 below. 


(Table T'.u 8) 




“ _"3 

-2 

— g 

""0 


" °b' 

^3' 

-t' • 

^1 

2 

Xg 

' ^4 


0 


29 

“*F” 

1 

0 

1 

2 


-2 

< 9 

2:o=-w 

2 5 

- 2 

5 

1 

0 

-1 

5 



^ ' iM ” 

*5 

■ ■ ' 19' ■ 
™ . 

0 

0 

.... ™ 

5 
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st ep 4 i Dual simplex algorithm gives the optimal solution 
as in tableau 9, 


(Tableau 9) 



c . 

0 

“I's 

. B 2' ' 

O"'’* 


■ Cb 

-sr 

■ X- 

^2 

-"*3 

^4 

0 

X2=21 . 

-5 

-2 

1 

0 

0 

x^=9 

2 

-5 

0 

1 


35 

3 

2 

0 

0 


Xte,rat.i pn 2 

St ep 2 s All the variables are within their bounds. 
Optimal solution is 2=35. 
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ASSIGNMENT POLYTOPE 


5.1 


Introduction 


The assignment problem of order n is 


n n 

min Z 2 c.- ^ x. - 

i=l 3=1 


subject to 


CAP) 


n 


2 X- . = 1 , i=l,2, ...jn. 

j=l 


n 

i=l 


X. . - 0 or 1 

Associated with this (LPP) Linear Programming Problem 
(without integer restrictions) is the assignment polytope 
(APj^) of order n which is the set of all feasible solutions 
to the assignment problem. Justification for the formulation 
of assignment problem as LPP lies in the fact that the 
extreme points of AP^, thought of as nxn matrices are 
precisely the permutation matrices (assignments). For 
proof refer to ' [ 92 ] • This also follows from the concept 
of total unimodularity of matrices [29,58]. Many methods 
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are available for solving the assignment problem. See 
for example [11,49,98,117 ] , 

The structure of the assignment poly tope helps to 
imderstand the assignment problem more clearly. A general 
study of polytopes is given in [ 71 ] , Klee and Witzgall [96] 
have discussed various aspects of transportation polytopes 
of which assignment polytope is a special case. However 
there are many results which are either not true for the 
general transportation polytopes or are not known. 

In this chapter jwe study the adjacency of the 
vertices of the assignment polytope, the structure of its • 
faces and some properties of the assignment polytope. 

Our approach is completely elementary and simpler than 
the work of some other research workers in the same 
direction [ 12 ] , 


5.2 


A d jace ncy „Re_s ultg 


Before we derive the adjacency results we need 
some definitions and notations . 

With any permutation ir £ we associate 
a permutation matrix xCtt) = (tt)^ as follows : 

Xj^^Cir) =1 if irCi) = j 


0 otherwise. 
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Thus, whereever convenient, vre may identify permutation 

matrix xCtt) with permutation ir only. Furthermore, we 

may use sometimes the notation ijig....!^ to mean the 

permutation / 1234 \ 

Siji2....i^J \2143./ 

either be written (in usual notation) as (12) (34) or 

as 2143, This notation in which the first row in its 

natural order is understood and hence not written is 

sometimes referred as "two rowed notation". Let 

TT and cr be any two permutations. Define 

I = ■^i ,0 ) |x^^ (tt) = 0 and ( a ) =0^ 

Let c - be an nxn matrix such that 

~ 1 if (ijj) 6 T 

= 0 if (i,j) ^ I 


Define the hyperplane H as 


n n 

H j S 2 c.. 

i=l 3=1 



= 0 


which may shortly be written as H : cx = 0 ( the 
manner of multiplication being understood). The line 
segment joining x(7r) and x( cr ) is denoted as |x(Tr),x(cr )~| 
or simply as [ v, 0 ] , 

Theorem 1 : H is a supporting hyperplane to AP^ 

containing x(Tr) and x(cr ), 
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Proof : For any x = ^ 

CX = S C. . X;. • > 0 

i,3 - 

and cx = 0 for x = xCtt), x(o ) (by definition of c). 

This completes tho proof. 

Using any tv;o permutation tt and cr , we 

can construct a supporting hypcrplane as above* 

This hyporplane will be used throughout this chapter 

(and tho next chapter also). The intersection of this 

supporting hyperplane (w.r.t. rr and e ) with AP^ 

i.G. HA AP is referred as the face AF(Tr,cr ) of AP„. 

n — n 

We observe that AF(tr, o ) is convex and tt and <y 
are in AF (tt, c ) . 

¥e now prove a theorem which characterises 
the adjacency of any two permutations on AP^ in terms 
of their adjacency on AFCv, cr ). 

Theorem 2 : Tuo permutations ir and d are adjacent on 

AP^ iff TT and cr are adjacent on AF (ir, a ), 

Proof : If TT and cr are adjacent on AP^^ then [■m, o] 

is an edge of AP^, 

Any point x 6 [tt, o] has a unique representation 
as a convex combination of the vertices of APj^ and that 
unique rc-prosontation is given by 

X = ocxCtt) + (l’-oc)x( a) for some >=< such that 


0 S ^ < 1 
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==^> X 6 [ IT, or 1 has -unique representation as convex 

combination of the vertices of AF(ir, cy )• 

==^ [it, o] is an edge of AF(7r,q ). 

=v> [ VjcT ] are adjacent on AFCitjct )• 

Converslv Let AFiir, a ) have m vertices^ namely, 
such that 1^= v, cs and rr and cr 

are adjacent on AF(v,o ), This implies any point x 
on the line segment Ivjcr] has a unique representation 
as a convex combination of the vertices of AFCvycy ) 
and that unique representation is given by 

X = pC^x(ir) + ^ ^ '^1- ^ 

* 4 “ ziz *1 ^ 

^ ^ nf • n| 

Let X = s' c(.x(T.-) , o(. > 0 , S 0(5^ = 1 

i=l ^ ~ ’ i=l "■ 

...( 1 ) 

bo any representation of x as convex combination of 

vertices of AP^ -where vortices of are named as 

TijTgj • • • jTni such that = v, Tg = a • Since 

X lies on the hyperpiano H, therefore 

n ! m n j 

0 = cx = S °<^cx(T^) = 2 <=(. cx(T. ) + s' cx(Tj|) 

i=l ^ ^ i=l ^ ^ i=m+l 

...( 2 ) 

Since H is a supporting hyperpiano to APj^ and T^ S' AF(v, a ) 
for 1 > m, 

cx(Tji^) > 0 for i = m+l,..,,n. 
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Thus in order that (2) may hold ~ ^ ~ Ki+l,,,,,n 

and this reduces (1) to the following form 

m m 

x= 2 x(T. ), °(±> 0, 2 = 1 

i=l ^ ^ - i=l 

NoWj since tt and a arc adjacent on AFCtt, cr )} therefore 

X = '=<T x(Tr) + °c„x(a ) (recalling that tt 5 Tn , 

^ ^ c E T^) 

i.c, X in (1) is expressible as a convex combination 
of x(t) and x(cy ) only which implies tt, cr are adjacent 
on 

The importance of the above theorem lies 
in the fact that adjacency of any two vertices on 
AP^ can be decided by deciding their adjacency on AFCtt, u ) 
which consists of fewer vertices of AP^^* The following 
theorem relates the nature of tt and cr to the number of 
vertices of AFCtt, cr )* 

Theorem 3 j If ir and a arc such that is a 

product of p disjoint cycles of length > 1 then 
AF(ir, cr ) contains 2^ permutations (including tr and 0 )« 
Proof : Let the p disjoint cycles of oir’i be 

^l? ^2’****^p 

such that each = l,2,...,p has length greater 


than one. 
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Let S = 'U Sj,. ‘ 

keQ 


where 


Q e 



i(i)ji e 

Construct a permutation T x^rith respect to set S as 
follows : 


Define J 


V" 


T(o) = cy (d), if d e J 
TCd) = yr(d), if d S' 1 

We show that T G AF(ir,(j )• 

I = |L,2,...,nJ. and let K = ^ G I U (d ) = ^(d >3^ 
i.e, set of indices corresponding to cycles of length 
one. Then T(j) = cf(d) = ir(d) j G K. 

Also for d ^ Jj have TCj) = cy(d) snd for 5GI~(JUK;)> 
we have T(j) = wCj) 


2 Sc.. x_. . (T) 

n T X J X J 


2 2 c • . . (T ) + 2 

IGF. d igj 


2 c. . X. .(T) 
1 10 


+ S 2 

iei-( jUk) j 




and each of the three summations on right hand side 
is zero separately. 

Therefore cz(T) = 0 which implies T G AFCtt, a). 

For different choices of set S we can construct 
different permutations T belonging to AF(7r, cr). Since there 
are 2^ choices of Q and each choice of Q gives a permutation 
T, we got 2^ permutations in AFCtt, cr ), 
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To prove that these are the only permutations 
6 AF(-f, ct)) let T he any permutation such that cz(T) = 0, 
Clearly cx(T) = 0 should hold separately for each row 
2 of matrix c and x(T) (by definition of c) 

=> T(;i) = cT(j) = -rrCo) for o £ K. 

and T(j) = crCj) or 7r(3) for j £* K. 


Let 

J = 

^1T(3 ) =(J (d)| 

Thus 

for 

3 G I- (JUK ) , T(3) = v(3) 

Let 

S = 

{a( 2 )\d e Jj. 

If 


|'ir(3)l3 G j'J- 

*- »»v/ 

Then 

S = 

S , because if not, then there exists 

1’ 


at least one j £ Jjsay, j = tj, such that S 

and cr (t^) 0 

=> T(t3_) e S and T(t^) ^ (because T(t^) = ff(t^)) 

=> T(t^) Tr(j) for 3 G J 

=> T(t^) = irCj) for some 3 ^ ) 

=> T(t, ) = TCj) for some j G I- (J4JK) (because T(j)=-ir(j) 

for 3 G I- (jyK); 

— > T is not a permutation. Hence S = S^. 

ilso 0(3) 7^ 11(3) for 3 G J. Therefore S is a union of 
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disjoint cycles of such that each cycle is of 

length greater than one i.e. S = U Therefore 

h£Q 

any T £ AFCtt, cr ) is obtainable from some set S of the 
above type only which implies AFCir, a ) has precisely 
2^ permutations and they can be constructed by the 
knowledge of the cycles of 

Bemark 1 i It can be observed that 2^ permutations 
of AFCtt, cJ ) can also be generated by premultiplying 
w by cycles of the combinations of cycles 

being taken in all possible ways. Thus members of 
AFCitjO ) are 

3 f 3 j 3 j 3~ 3^7T j 3<i^ 3^7r ^ ^ S*. 3 tt ^ 3^3 tf 

and these are precisely 2^ in number. 

Corollary 1 ; An immediate consequence of the above 

Theorem is that if ir and cr be two optimal assignments 
such that cr-jr”^ is product of p disjoint cycles of 
length greater than one, then there are at least 2^ 
optimal assignments. 

The following Theorems (4,5 and 6) 
decide the adjacency of two permutations and the 
number of vertices (order of adjacency) adjacent to 
a given vertex on AP„. 
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Theorem 4 j If tt and cr he any two permutations such 
that ott"! is product of p disjoint cycles then 2^ 
permutations £ AF (it, o ) can he so paired such that 
2^"^ line segments obtained hy joining these pairs have 
the same mid point as the line segment [x(Tr)jx(c 5 ' )] • 
Proof s Let G AFCir, e ) he obtained through 
Q C ^1,2, ...,pj^ as in Theorem 3. Let Tg £ AF(ir,cr ) 
he similarly obtained through Q ’ = 

Then obviously j 

I x(T^) + I xCTg) = I xCtt) + I x(ff) 

This completes the proof. 

Each of the 2^^ pairs of assignments obtained 
in Theorem 4 are referred to as complimentary pair of 
assignments ( or simply complimentary assignments ) w.r.t. 
the given pair of assignments v and c • 

Theorem 5 : Two permutation tt and a are adjacent 
on AP^ iff cfTr“^ is cyclic. 

Proof : We know by Theorem 2 that two permutations 
TT and o are adjacent on AP^^ iff they are adjacent 
on AFC-ir, cr ). This is so iff AFCir, p ) contains just 
two permutations (tt and <7 itself) which happens iff 
e TT*"^ is cyclic because otherwise ir and b are 


non-ad jacent in viev; of Theorem 4. 
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Theorem 6 : The nmher of extreme points adjacent to 

any given extreme point on AP^ is the same and equals 

n 

N(n) = 2 (r-l) ! n 

r=2 ^r 

Proof : That the order of adjacency for each vertex 

is the same is obvious. To compute this number, 
let TT = 1,2, , . . , r, , . . ,n be any permutation in two 
rowed notation. If cf is adjacent to ir then a must 
be of the form so as to form a cycle of length r with 
IT (o written in two rowed notation). There are 

choices of cycles of length r. Keeping one element 
fixed, any such cycle can be permutated among its 
olemonts in (r-Dj ways to give rise to different 
permutations. Also since minimum length of a cycle to 
give rise to cr -r is 2, therefore 

n 

N(n) = 2 (r-l)! n^ 

r=2 ^ 

Here below we give some numerical estimates of the 
adjacency which, as clear from the tableau, is of 
very high order even for the assignment polytopes 
of small order. The table also gives the dimension of 
APj^. We know that nxn assignment problem when written 
as linear program has 2n rows out of which any 2n~l rows 
are linearly independent [ 76, Chapter 9, Transportation 
Problems ] and since APj^ lies in space, is linear 



122 


vector space of dimension n^“(2n-l) = (n-1)^. The 
following table gives the value of dimension of 

AP and N (n ) for 1 < n < 10 . 

n - - 



Dimens ion 

No of ext.pts. 

No* of adjacent 
ext.iDts. 

1 

0 

1 

0 

2 

1 

2 

1 

3 

4 

6 

5 

4 

9 

24 

20 

5 

16 

. 120 

84 

6 

25 

1 720 

409 

7 

36 . 

1 5040 

2540 

8 

49 

40320 

16064 

9 

64 

362880 

125664 

10 

81 

3628800 

1112073 


(Some Estimates for AP^) 


In the following Theorems,? and 8, we mention some more 
facts about the face AF(ir, o ) which helps us to imderstand 
the s tructure of more cloarly. 

Theorem 7 : If each cydo of av-"! is of length, k, then 

the distance botwoon any two extromo points adjacent over 

1/2 

AFCvy a ) is (2k) , 
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Proof : If and be any two adjacent permutations 

over AFCtt, cy )j then T is cyclic. Since T2_ and T 

J- ^ 2 

arc obtainable from the Imowledge of cycles of air 
therefore 12.^2^ is some cycle of atr"^ whose length 
is given to be k, 2 . 

Obviously llx(T2_)-z(T2) |1 = ( 2 k) 

Theorem 8 s If erv”^ is product of p disjoint 
cycles then any permutation £ AFCtt, a ) has exactly p 
permutations 8 AF(v, a) which are adjacent over AFCv, c )• 
Proof : Let and T2 be two permutations over AF(ir, cr). 
They are adjacent over AFCtt, c ) iff ^^Tg^ is cyclic. 

Also cycle of 12.^2 cycle of o v , Obviously 

such choices of T^in AF(v, a ) for a fixed T^ are 
precisely p. Hence the proof. 

Geometrically vjo can think of face AFCir, a ) 
as in Figure 1 below, Tho figure also demonstrates 
the complimento.ry character of T and T such that 
mid point of [ x(T2_),x(Tg) ] being tho same as tho 
mid point of [x(v),x(a )J. . 
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However all the faces of are certainly not of 
the kind considered above. For example j there 
can be faces of with only 3 vertices on them 

as the following example shows. 

Example 

jr = 123456 

<7 = 134526 (Two rowed notation) 

T = 134256 

¥e observe Tra”^ = (2345), cr T~^ = (52), irT”^ = (234) 
so that ir, cr and T are adjacent to one another over 
AP^ giving rise to a triangular face of AP^ (Figure 2) 

T 



(Figure 2) 

¥g now give below few properties of the 
inverses of the permutations and the relation between 
their corresponding faces. 

Clearly a permutation ir and its inverse 
7r“^ are adjacent. 
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Theorem 9 : If and o be two permutations 

adjacent over AP^ then so are tt""^ and cr~l. Also 
distance between and is the same as distance 

between tt and a . 

Proof s Let ir •, irCl) ,Tr(2), , ,Tr(k) , .jirCn) 

e i O’ (1 ) j d(2 ) j » • • jc? Ck ) } • • • j cj Cn ) 

Let I = |i : rCi) o (i)^ 

Without loss of generality) let I = -^,2,,..)k^ 

so that 7r(i) = o (i) for i = k+l,..,,n. Again, we 
assume without loss of generality 

o C 1 ) ” m ( 3.*^ 1) , 1 —" 1)««») k^ 1 

o (k) = -n-(l) 

For i I, o (i) = Tr“l7r(i) = i 

For i G I, Tfl o(i) = Tr"^ir(i+1) 

(i = 1, • . . ,k-l) 

= i+1 

and o (k)= 1 

s=s> vT^a is cyclic of length k* Thus tt"*! and cr~l 
are adjacent over AP^. 

Obviously jlx(Tr”^)-s( o || = | jx(Tr)-x(o ) { 1 

Corollary 2 : Permutations ir and o are adjacent 

over AP iff and cr”! are adjacent over AP„* 

Proof : If V and. o are adjacent, then each of 

o 'nr~l,'r o"*^, 7r'“^cT j o is cyclic. If tr and o 
are non- adjacent then Tr"^ and o“^ are also non- adjacent 
aid accordingly each of cyv’'^,Tr cr*’^,v'‘^cr , cr“^ is non-cyclic 
and is product of same number of disjoint cycles. 
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The proof now follows obviously. 

The following theorem establishes the 
relation between the faces corresponding to inverse 
permutations • 

Theorem 10 i If tt and cr be two permutations 
over AP^ and AFCtt, ct ) and AFCtt""!, a~I) be faces of 
APn corresponding to permutation pairs ) and 

( 7 r“l, cy”^) respectively, then 

(a) IafCtt, cr )[ = IaFCtt"!, a“l)l 

(b) [AFCtt, cy ) ]“^ = AF(7r“l, cfl) 

where [AFCtt, cr ) ] "^ denotes the set of inverse 
of permutations in AFCyf*^, )• 

Proof I Part (a) is obvious because if ■jr'“l is 
product of p disjoint cycles then so is c-Itj-, 

If IT and a are two adjacent permutations 
then part (b) also follows obviously. So we assume 
that IT and cr are two non-adjacent permutations so as 
cr is product of, say , p disjoint cycles of length 
greater than one. Let T^ G Aj.' (yr, o' ) be obtainable 
through S such that 

® = l^Q s {1.2 p} 

J = ^yr-^(i)li e sj 
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Let S’ = I-S and J* = 

Obviously jf\ j* = 9 and J U = I 
Thus T^(j) = a(o) ,361 
T^(o) = Tr(j) j 0 G J* 

Let = -^* e Jlv'‘l(j)ejJ and Jg = -(jeJl ^*“1(0 ) ^ 

We observe J_ f\ = 9 and J, t = J 

r ^ ^ ^ tr A 

Let Jg = -^‘GJ*l'ir“^(j)GJ‘*|" and = ^GJ* 

We observe J 3 A J 4 = 9 and J 2 \_j ^ 

Now define T^GAFCtt"!, cj"^) as follows s 

T^Cj) = 0-1(3) , 3eJ-3_UJ4 

t*( 3) = v-i(3) joejgU J 3 

^ *•* 1 

ftTe prove in the following that T^ = Tj_ 

For 3 GJ^,T^T*( 3 ) = T^Ca'lCj)) (for 3eJ^3^“^(o 

= a (^”^(3)) => a“l(3)GJi) 

= 3 

For = T^(7r“l(3)) (for 3GJ2,Tr“l(3 )^J 

= tt(v “1 ( 3 )) — > 'rr"l( 3 )GJ*) 

= j 

For 3 GJg , T jT* ( 3 * ) = T^(Tr^l( 3 )) (for 3 GJ 3 ,Tr“l (3 )GJ*) 

= v( 7 r“l ( 3 )) 

= 3 

For 3GJ^,TT*(3) = T( o“l(3)) (for 3eJ4,TT^^(^ )^J* 

= O ( CJ-1(3)) = ^"^(3)0/ 

= j = <’ ‘^(Uej ) 
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which proves = t“^. 

Since lAF(ir“^j cr"l)j = {aFCvsO )ij (by part a) 
the proof of part (b) is complete, 

5.3 Po lytope _ 

The ^^ifft anc6_^ ji (^aj^ b_) between a pair of vertices 
a and b of a convex poly tope is the minimum number of 
steps needed to reach vertex b starting from vertex a 
and always moving along an edge of the polytope, (A step 
means moving from one extreme pojnt to another adjacent 
extreme point along an edge of the polytope). The 

d of a polytope X is defined as the greatest distance 
between any pair of vertices of the polytope X i.e, 

d = max dCajb) 
ajbCX 

Theorem 11 : The diameter of AP is two, 

■ ■ n 

Pj’.osl'. ’ To prove that the diameter of AP^^ is two, we prove 
that if TT and a be two non** adjacent vertices of AP^^ it 
is possible to cons tract a vertex T of AP^ adjacent to 
both TT and o . 

Since v and a are non- adjacent, is 

product of, say, p disjoint cycles S^,S , ,,.,S . 

If we constract a permutation T as in Theorem 3 
by taking S to be one of the cycles S^jS , ,,,,S , then 

j- 2 p 
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clearly irT"^ aiid arc cyclic which implies T is 

adjacent to both v and a , In general jwe can get 2p 
such permutations (for p^ 3) aujacent to both tt and ff* • 
IIo^^reverJ this number is not enhaustive. 

In [12,126 ] 5 the diameter of AP^ is showi 
to be equal to two. In [l2l] jthe diaoieter of the assignment 
polytope, traveling salesman poljT-tope end polytopes associated 
X'Tith a large number of combinatorial problems has been 
sho'vm to be less than or equal to two. In the case of 
assignment poly tope we have achieved the similar result 
above in a straightforward manner* 

5.4 Graph of G(APj^) 

Before we deal with the main result of this 
section, we give some definitions [78] of graph theory. 

The paths in a Graph G are called edge~ disjoint 
paths if they do not have any edge in common. A graph 
Map. -N-conne cted provided there exist N mutually edge- 
disjoint paths between every pair of nodes of the graph G. 

The distjmee d(a,b) between a pair of nodes a and b is the 
length of shortest path containing them. The diameter D(G) 
of a connected graph G is the greatest distance between 
any pair of nodes i.e. 


D(G) = max d(a,b) 
a, b6G 
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B s ki ^ s C on .1 c t ur e [ 12 ] 

If G(i\P^) is this graph v/hose nodes are the. 
vertices of and edges are the lines join-ing the 
two adjacent vertices of AP^, then G(AP^) is N (n)~ connected 
where N(n) given by 

n . 

N(n) = 2 n_ (r-Di 

r=2 r 

is the number of vertices adjacent to any vertex x on 
APn * 

of the Conjec ture' [ 89 ] 

G(AP^) is 8dge-l'I(n)-connected, 

We will first prove this weaiker form of the 
conjecture and then a small observation will extend the 
result to all polytopes of diameter twoj a large class 
of which is considered in [ 121] . 

Proof i For any vertex a "of. APj^j let 

X(a) = G APj^ ; X is adjacent to aj* 

For any two non adjacent vertices a and b of AP let 

n 

X(a)n X(b) = KCajb) 


’'•‘Presented at Ninth Annual Convention of Opera.tion Research 
Society of India held at Calcutta (India) in Dec, 1976 and 
was awarded the prize for the best student paper [89]., 



131 


K^(a.jb) — iC(a)“K(3.j'b) 

K^(a,b) = X(b)-K(a,b) 

There are three types of paths as described below . 



IC ( a j b ) 

(Figure 3 ) 

(I) For k G K(a,b) 

PCajkjb) = [ (ajk) j (k,b)] ...(3) 

is a peth which joins a to b. 

For k^G IvCajb) and kg G K(a,b) 
the patt-B P(a,k 2 _,b) and PCajkgjb) are two disjoint paths joining 


a to b. 
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Thus 


IS a 


= [PCajkjh) J k £ K(a,b) ] 
set of disjoint paths joining a to h. 


Let 


= [(tjjtg) 

0 -h 

£ (aj b ) } 

tg e Kg 

(a,b). 



h 

is adjacent 

to tg] 

« 

s* <c. 

S S 

iuch 

that 




>li 

e s 

and 

(t3_,t^) £ S* 

=> tg 

= *2 


e s 

and 

(t;[,tg) £ S* 

=> tj_ 

= H 


For ^ ^ 


PCajt^jtgjh) = [(a,t^), (t^jtg), (tgjb) ] ...(4) 


is a path which joins a to h. 

For £ S*, (tg.t^) £ S*, the paths P(a,tj^,tg,b) 

and P(a,t 3 ,t^,h) are disjoint paths joining a to b. 

Thus 

Sg = [PCajt^jtgjh) : ^ ^ 

is a set of disjoint paths joinuig a to h. Observe that 
each path of Sg is also disjoint from each path of S^. 

Define 

S* = [t^ ! (t^.tg) e S* lor some t^ £ KgCa.h) 3 
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Let 

K = K^(a5b)-3^ 

P = IC 2 (a 3 b)~S 2 

Obviously 3 K and P have the same nuraber ox elf-ments and 
for any 6 Kj there exist no ^2 G P such that 

t2_ and t^ are adjacent* 

r (5 ) 

III Now for T^’e K, chooses any t^ S P 
Since ■b2_j’t2 are non adjacent 

tit”^ is not cyclic. 

Now every permutation which is noo cyclic can 
be written as a product of two cucles [ 12 ] i.e. if 
is product of k disjoint cycles 


TT 


(v(l) 5 . . * ,v(n. )) (irCn-i+l) > • • • 


then 


TT = (ttCI) , . ,v(n) ) (v(n) . • •’^[^2^ vCn^^)) 


= uv 


V7here u = C 7 r(l)...v(n)) and v= (irin) , ^ 
are both cyclic. 

Let 


1--1 
a ‘^3 


tif^-" = U3_v^ 


where and v^ aro cyclic 

Obviously is adjacent to both t^ and t^. 


.»v(n)) 
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Then 

is a path joining a to b . 

We observe that 

v^^tg If K U P 

(a) If v^tg e (K^(a,b) IjK^Cajh))', then we 

(b) If v^t^ e (K^(a,b)-K), then path (6) is 
from path of the type (4), namely, 

(ajV^tg), (Vjt^jt-i^), Ct-] 2 >b) 



where t^^ is adjacent to £ (K^Cajb)-?) and t-]^ t^ 

otheriwse ^ ^ ^2 

r36=)> tg ^ P 

and this contradicts the choice of tg made in (5). 


(vj^tgjtg), (tg,b) J 
....( 6 ) 


are through 
dis j oint 
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(c) If 0 (K^(a 3 b)--P )5 the argument is similar 

as in case, of (b) and (6) is a path disjoint to the 
paths of the tj^'pe (4), 

For tg 0 K, t^ 0 P, the paths 

PCajtjLjV^tg, tgjb) and P(a, t^jV^t^j t_^jb) 

are two disjoint paths joining a to b. 

Thus set Sg = [PCajt^jVjtgjt^jb) s t^ 0 K, tg £ P 3 
is a set of disjoint paths joining a to b. Observe 
that each path in Sg is disjoint, from each path in 
Sj and S and the number of elements in S 0j S S 
e quals N (n ) , 

If a and b are adjacent vertices than (a,b) 
constitutes one path and the procedure outlined above for 
non-adjacent case can be repeated by replacing X(a) 
and X(b) by X(a)- -^b}- and Xih)- respectively. 

E xampl es 


For the purpose of demons trati on, the following 
simple examples of the assignment polytope are considered. 

For n < 3 we have paths of length one only. 

For n = 4, it can bo easily seen that any two non-adjacent 
vertices tt and cf on P^ are such that 


XCtt) = X(a ) = K(v,a ) 
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• li'i “the above proof, the number of vertices 

adjacent to any vertex on the assignment polytope AP^ 

is the same. However, this may not be true in general. 

case of any poly tope P , 

Define N = min N(a) 

aep 

Then G(P) is edge- N- connected. The proof will be 
exactly similar. 


In 



CKAPTER - 6 


TRAVELING SALESMAJ M POLYTQPE 


6,1 Introdu ction 


The traveling salesman problem of order n is 
n n 

min 2 2 ^ x. . 

i=l j=l 


(TSP) 


subject to 


n 


2 

j=l 



1 , i 


1,2 


> • • j 


n 


n 

2 

i=l 



= 1 » 


j - 1,2, 


# • ^ ri 


X. , = 0 or 1 
13 


and the solution set of those x> .»s for which x. . = 1 

^0 13 

forms a tour where a tour is defined as the cyclic 
permutation (i. ,i , ...ji ) such that the solution 

»-*- o Jl 

X. . = X,. ^ = ... = X,- ■? =1 and zeroes elsewhere 

^ 1^2 ^2 3 nn 

and vice versa. Evidently there are (n“l)| tours which 
are all feasible solutions of (TSP). 

Associated with this traveling salesman 
problem of order n is the traveling salesman polytope (TP^^) 
of order n which is defined as the convex hull of permutation 
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matrices corresponding to traveling salesman toiirs 
(or tours simply). Evidently TP^ is a collapsed 
polytope of AP^* Murty CllS] has discussed the 
adjacency properties of extreme points of collapsed 
polytopes. If two extreme points are adjacent over a 
polytope, they are also adjacent over its collapsed 
polytope. However the converse is not true. Thus if 
TT and a are adjacent over AP^, then tt and a are 
adjacent over TP^. 

Many exact as well as heuristic methods 
for traveling salesman problem have been developed. 

See for example [39,49,58,79,80 ]. However, 
in all the existing methods, the computational time 
increases ex p onen t i ally with increase in the size of 
problem. In fact it has been shown that the traveling 
salesman problem is EP^h ^d, . Refer to Karp [ 93 ] 
for details. However, Rosenkrantz, Stearns and Lewis II [132] 
have considered some polynomial time algorithms for 
finding ” i: approximate optimal tours", but not 
necessary optimal tours, for the traveling salesman 
problem. 

One good reason that good methods could 
not be available C£an be that very little is known to-*day 
about the traveling salesman polytope. In this chapter, 
we make an attempt to study the traveling salesman 
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polytope* Primarilyj-we conc-ern ourselves with the 
adjacency of the tours on approach developed 

for AP^ in chapter 5, is further explored over TP^ to 
know the adjacency of the tours on TP^* study, 

however, is incomplete* 

In section 6.2, we derive some adjacency 
results over Section 6.3 concerns with some 

adjacency rules and their applications for the development 
of heuristics. Section 6,4 lists some problems for 
future research, 

6*2 Adjacency Results 

Corresponding to tours tt, p we define the face 
TFCtt, p ) of TPj^, similar to face AFCtt, p ) of AP^^ as follows : 

TFCtt, p ) = ; cx(T) = 0 and T is a tour^ 

where the matrix c is constructed w.r.t, tour ir, p 
as in chapter 5. 

For a tour tt, let tour o be such that p tt”! 
is product of p cycles each of length > 1, Let 
S = Sj_US2 ***USp and let jSl = k, then we say 
p G simply p G Z^. Note that Zg = 9. 

In view of the folloxiring Theorem (similar to 
Theorem 2 of chapter 6) the adjacency of any two vertices 
of TPj^ can be decided by deciding their adjacency over 
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TFCtt, a ) which contains much smaller number of vertices 

than TP , 
n 

The ore m 1 : Two tours it and a are adjacent over TP^ 

iff they are adjacent over TF (ir, )• 

Proof : Exactly similar to Theorem 2 of chapter 5, 

Theorem 2 i If [ TFCir, a ) 1 < 3, then tt and o are 
adjacent over TP^. 

Proof i If TFCtt, a) ~ 2 i,e» TFCttj o' ) = 
then obviously [■?r,o' ] is an edge of TP^. 

Let TFCtTj c 7 ) = -(TjTt, (7 V • Since T,Tr and tj 
are three vertices of TP^^ (lying over a hyperplane 
H t cx = 0) they cannot lie over a line. Therefore 
[ TTj cr ] is an edge of TF(7rja ) which implies it is an 
edge of TPj^. Hence the proof. 

In fact we have also proved above that [TjTt ] 
are [ T50] is also an edge of TF(vj o ) so that all the 
three tours of -fT,v,cr'\ are adjacent to one another, 

'v J 

It may be remarked that the above theorem 
indiccites the existence of triangular faces of IP^* 

The following example shows the existence of such 
triangular fac^s on TP^, 

Example 1 

Let tr = (1 2 3 4 5 6 7 8 9 10) 

0= (174936258 10) 
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Then TFCr, cr ) = i'lr, cy,T^ 

where T= (123674589 10) 


T 



(Figure 1) 

Th eorem 3 : If for tours tt and Oj is product 

of two cycles such t hat at least one of them is of length 

two then ir and are adjacent over 

p roof i Since is product of 2 cycles j 

lAF(Tr,o )1 = 4 

while TT, a are two permutations in AF('r, a )? other 

two can be obtained through set only; which consists 
of elements of cycle of length 2« (because a permutation 
obtained through S can also be obtained through S* by 
defining T^(j) - Tr(j), j ( J» and 3 2^ J* 

where J’ = /v“*^(i)li £ S» V . Since = 9 , these two 
permutations do not lead to a tour. Therefore TF('Tr, o ) 
contains ir and cr only which implies ir, o are adjacent 


on TP . 
n 



143 


C orollary 1 : A3L1 tours on TPg are adjacent to one 

another. 

Proof : If TT and cr are permutations on 5 symbols then 
a if not cyclic, can have either two cycles of 

length two or one cycle of length three and one cycle 
of length two. In each case, the proof follows by the 
above Theorem. 

The Theorem 2 can be interpreted as saying that a 
necessary condition for tt, o to be non~adjacent and TP^ 
is that ITFCtt, cr ) | > 4. 

However the above condition is not sufficient 
as illustrated by the following counter example. 

Exampl e 2 t Let tt - (1 2345678 9), 

cr = (1 47 639825), Using two rowed notation, the 

elements of the face TF('Tr, cr ) are written as below : 

base row 1 2 3456789 
v! 234567891 
a i 459713628 
cr = (2 4 9) (1 8 6) (3 5 7) 

Elements of TF(7r, cr2_ (other than rr and ) are s 

T^j439567821 = (1 4567823 9) 

T„:25476389l = (125634789) 

1^:459763821 = (1 4782563 9) 

3 **<!*.•** 
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Now by observing tt, a ,T ,T ,T <5 and their 
permutation matrices it is oberved (because of underlined 
positions i.e, any one of the rows 5th, 7th or 9th of the 
corresponding permutation matrices) that line segment- 
[ x(ir) ,x( cr ) ] has a unique representation as a convex 
combination of extreme points of Thus even tho\igh 

1 TFCtt, cr ) 1 5, TT and a are adjacent. 

Geometrically the face TF('7r, a ) can be viewed 



Complemen tary to xrrs : Two tours T^ and Tg on face 

TFCtt, cr ) are called ” complementary tours” w,r,t, tt and cr 
or simply ’’complementary” if they can be obtained 
through sets S and S’ such that S - (J S, , S’ = I j S, 

, ^ keq> 

where Q Cl -iljS, and p is the number of disjoint 
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cycles of a of length greater than one. Thus if 
Tj_ and are complementary w,r*t* ir and a then 

1 xCtt) + i 2 J:(cr)= 2 ^ x(T. ) + i x(Tp)« 

2 2 2 1 2 “=^ 

¥e now state a _suffj._c_ ient condition for tours 
TT^cr of TP^ to he non-adjacent . (on 

S uffic ie nt condition : Two tours tt and o are 

non-adjacent on TP^ if|TF(7rj o )1 > 4 and TF(ir,c)i ) contains a 

pair of complementary tours* 

Proof : If TFCitjo ) contains a pair of complementary 

tours then ir and o are obviously non-adjacent and 
iTFC-rr, ct)1 > 4. 

Howeverj the above condition is not necessary 
i,e* Sven if TFCtt, a ) does not contain a pair of 
complementary tours, ir and a may still be non- adjacent. 

This is shovjn by the follovring counter example. 

Let TT = (1 2 3 4 5 6 7 8 9), a =(1 5 9 4 8 3 7 2 6) 

Using two rowed notation, the elements of face TF(ir, p ) 
are written as below : 
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base row ; 123456789 

■nr : 234567891 

<y : 56789123 4 

=(2 5 8) (3 6 9) (1 4 7) 

.Elements _of TF(Tr, a ) (other than tt and o ) are : 

T^: 534867291 =(156723489) 

T : 264597831 =(1 2678345 9) 

2 

T i 237561894 =(1 2378945 6) 

3 

None of the pairs (T 2 _,T 2 ), and (TgjT^) is complementary 

w.r.t. IT j cr • However j it can be observed that 

I xCtt) + I x( a ) =1 x(T 3^) + I xCTg) + 

which establishes our claim. 

Geometrically the face TF(ir, cr ) lool?S JLitee 
the following Figure 3, 


T 

3 



(Figure 3) 


CO 5m 
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Ihe follo^^ring theorem gives necessary and 
sufficient condition for non-adjacency of two tours on 
TP^ in the special case when [TF(Tr, a )| = 4, 

Theor em 5 : If|TF(irjcr )(= 4, then tt and a are non- 

ad Jacent on TP^ iff TFCttjCT ) contains a pair of 
complementary tours. 

P roof ; If lTF(ir, c?)! =4 and TF(ir, or ) contains a 
pair of complementary toirrs, then tt and a are obviously 
non- adjacent on 

, o ,T^,T 

be non-ad Jacent on TP^. If and Tg are not complementary 
tours 5 they can be obtained through two cycles S-j^ and 
such that 


Converse_ly let TF (tt, a ) 


(jr 




\ 


and TT.a 




u \ 




ke'Qg 


where and are subsets of P = -^^^,2, ,.,,p 
assuming that orv""^ consists of p disjoint cycles of 
length greater than one, such that one of the following 
holds 


(i) 

(ii) 


Q^U Qg is a proper subset of h,s,—,pj- 

We take each of the above cases separately 

Isi P* = P** ( Q2_ Li Qg ) ) let S* = pgp* ^ 


Define J* 
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Then T- and T„ ai’e iaich that 
1 2 

T^(0) = 7r(j) , D G J* 

TgCj) = 7r(j) , j £• J* 

ca(r) + a-c<)x(cy ) 7^ B2 c(T^) + a-R)z(T 2 ) for 0 < 

0 < B 

=> IT, a are adjacent which is a contradiction, 

C.ase (ii) Let S 

' ' pGQ^n Qg P 

Define = { Tr“l(i)li £ S* 

Then T^ and Tg are such that 

T^Cj ) = a (j ) , 3 £ J* 

Tg(j ) = c (j) r 3 £ J* 

==> oCx(Tr) + (l-°<)x(cy ) ^ Rx(T ) + (l*-B)x(T ) for 0 < ^ < 1, 

o *** 

0 < B < 1. 

=> TTja are adjacent which is a contraction. 

Thus Tj^ and Tg are complementary tours. Hence the proof, 
Bao [l26] has also studied the properties of 
traveling salesman tours by using graph- theory approach 
and has achieved similar results. 


Vj V| 



6.3 


Ir. .Aaei-J-Sa^ons 


6,3,1, Formulation o f Ad ;iac ency Rules 


Ue formulate below some rules to generate 
adjacent tour to a given tour, Ploweverj many adja-cency 
and non-adjacency rules can be formed over TP^ but it 
is difficult to say in any way that they are exhaustive. 

We are giving a brief account of a few of them only which 
are elementary in nature and have some applications in the 
development of heuristics for solving traveling salesman 
problem. Most of these adjacency rules are generated 
by permuting a subset of the elements of a tour among 
themselves by a suitable rule, 

iga is* .ip/i^_^^. . . ir+|5;+i 

be a tour. If k consecutive elements ir+l’^r+P* * '^r+k 
as marked are chosen and replaced by the sequence 
ip+gj ♦ * where 1 < s < k then the 


resulting tour is called a k- cyclic c hange of t. e,g. 
if V = (1 2 3 "?'5 6 7 8 ) then 2 3 4 5 when replaced 
by 345 2, 4523 and 5234 gives = (1 3 4 5 2 6 7 8) 

cJg = (1 4 5 2 3 6 7 8) and = (1 5 2 3 4 6 7 8) 

respectively. In this case each ” 1|2,3) is referred 

to as 4-cyclic change of ir. Geometrically such changes 
can be viewed over a cycle (Figure 4) in a clockwise 
direction. 
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(Figure 4) 

If the number of elements (i.e.k) is immaterial then 
the changes are simply referred to as cyclic changes. 

Rul6_ _1 : Any cyclic change of a given tour tt is a tour 

adjacent to ir on 

Prjx^f ; Without loss j of generality, let 

Let o be a k—cyclic change of rr i.e, 

d ” (l,^,.«.i, X"^s « • # i'^S"i, I'^k’^l , » • • , n ^ 

Define K = /kimCk) ^ a (k) 1, 

1. J 

Obviously K = -^iji+s-lji+k 
If I = then we observe that 

T(ir(i)) = cf (i) 

= ir(i) 



i 0 K 
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'^.■(•ir(i)) =0(i) 

i £■ K 

1 ^ Tr(i) 

Since a permutation ovor thrive bywbols ruicli that all 
the three syiiibols are mapped separatelvjis cyclic, 
therefore T is cyclic which implic?^ e is adjacent 
to TT. 

Vie observe that all such k-cyclic changes belong 

to Z • A special case of these cyclic changes when 

any two caisecutive elements of ir are changed is used 

by us to prove that TP^ is hauiltonian, (Section 6.3,2). 

Rule 2 ; If a tour a is obtained from tour ir by 

interchanging any two non- consecutive elements then cr 

is adjacent to rr on TP • 

n 

RJbQof • Without loss of generality, let 

TT = (1,2, . . j-1, j, j+1, ,, ,,k-l,k,Is;+l, « * ,yn) 

Let any two non- consecutive elements, say, j and k be 
interchanged ( 1 and n are also consecutive elements), 
so as to get 

<y = (1 j 2 j • • • , j~l, k, j+l, . • * , k-l, j , k+lj , * , ,n ) 

Let K = 'filTrd) a (i)_^*- 

Obviously K= / j-1, j ,k-l,k|- 

If T = then wo observe that 


T(Tr(i)) = a(i) 
= ir (i) 


for i K 
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cjid T(7r(i)) = cr (i) 

Tr(i) for i t K 

iJow T(Tr(j-l)) = <j (j-l) 

=> T(o‘) = k 

T(7r(k~l)) = <? (k-1) 

TCk) = 0 

T(7r(3)) = ct(3) 

==> T(j+1) = k+1 
TCirCk)) = (y(k) 

T(k+1) = 3+1 

I'iius T = is product of two disjoint cycles namely 

(3,k) and (j+ljk+l)* Now ty Theoi^eia 3, tt, a are 
adjacent* 


Tour cr as obtained above £ Z^» 

: If in a tour tt anj^ four consecutive elements 

are porniuted among theniselvos tbon the resulting tour 
is adjacent to rr on TP^^. 


EI'P-PX • ¥ithout lovss of generality, let 


Let 


wboro 


of tho 
Let 


TT = (1,3,3, ...Icjjc+Ijk+Pjk+S, ...,n) 
o — {l,3,.i>,»«»s, s+ 1 , s+2 , s+3, • « « , n) 

(s, ;'.+!, s+2,s+3) is some non- identity permutation 
symbols (k, k+1, k+2 , k+3 ) . 

J - ■^ilv(i) cy(i)} 
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liicn J call consist of at the most following five symbols 
i.e* J = *; k-l}kjk+l}k+2jk+3_y- 


If ± ~ a then 

T(Tr(i)) = 

oCi) 



Tr(i) , i ^ J 

and 

T(7r(i)) = 

cs (i) 



irCi) f i £ J 


This implies T(TF(i)) ir(i) for at most five values of 

i namely i = k-l,k,k+l,k+2,k+3* 

Thus T can be a permutation of the following types only 

(i) cyclic 

(ii) product of two cycles such that at least one of 
them is of length two. 

In ench case cr is adjacent to tt (Theorem 3). 

Since what is true of tour t = (lj2,...,n) 
is also true of any other torn.’ (simply the symbols are 
to be renamed) so for the sake of simplicity and 
convenience j in the following, we state results vrith 
respect to vs (1,2, ...jn)* 

Inve rsion t • ■ 

Let -ir .= (1,2, ...,k~l,k,T4a7i?^^T7r,l?^,’F^37PiP+l» a tour. 

Let le =. ( 1 , 2 , ,.,,k«l,p,p-l, p~2,..,,k+2, k+]jk, P+1,... jn) be 
obtained from t such that the block (k,k+l,k+2, . .. ,p-2,p«l,p) 
of elements appears in o in reverse order. Then ^ is 
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said to be obtained from ir by the inversion I(k,p) or 

simply a is Inversion of ir. 

The definition of inversion is due to Croes [5 9] . 

s Inversion of a tone tt is a tour adjacent to 

TT on IT . 
n 

P roof : If T = OTr“^ then permutation T looks like 

/ k ici-l k:+2 p-2 ,p-l p p+l\ 

yp.. p+1 k . k+l -**•., p-4 p-3- .p-2 P“l/ 

Wow there are two cases (i) p-Ic+1 = even, (ii) p-k+1 = odd 
(i) D-k+l = even : In this case T is cyclic and accordingly 
is adjacent to tt on 

Cii) D-k4-l = odd : In this case T is product of two cycles, 
namely 

[p ) p— 2 , p— 4 , • * • g , k*^^ , k“^2, k ) and (p"^l , p— 1 , p— 3 , « g g g , k’^3 , k*^! ) 

AF (it,'? ) contains four permutations two of which are 
v and cr , The other two (written in two rowed notation) 
are s 

base row ; 12 3 ..g. k-lgk ai:+iac+2.« .>.gP-a^.Pr2«nrl^^ n 

=2 3 4 .... p , k+1, k ,k+3, .,,P“4,p-l,p-2,p+l,p+2, ..1 
ig ~2 3 4 gagg k, p^"! j k"t'2 , k*!"! j g g , p— 2 , p— 3 , p , p— 1 , p"iT , g g 1 

It can be easily seen that both of and are non- tours# 

Therefore TFCitjO ) = 2 which implies ir and a are adjacent 

on TP « 
n 
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it ^lG 5 # Isi/ TT — I v^ 1:C*^1 t ♦ « • ^ l£"^s y « • 3!* « » D y p*i~ 1 • ^ p"i~ tj j • • ^ n } 

Let Ck,k+lj .. .jk+s ) and (p jp+1, . . .jp+t) be blocks of two 
eloiaonts of w separated by at least one element* Let a 
be obtained from tt by mutual interchange of these two 
blocks i.e* 

e ^ P"^l j j P"^t y ##*jrj**jk| k*^! ^ ^ k*^s 5 p*^t"^l ^ * ^n ) 

Then cr is adjacent to r on TP^. 

Proo f : For the sake of simplicity, we first assume 

that 1 < k and p+t<n. However, at the end, we will observe 
that such restrictions are not necessary. Let 


K = |il7r(i) a (!)}■ 


Clearly K = 

■{k-l,k+s,p-.l,p+tj^ 


If 

T = 0 7r“l, 

then 




T(Tr(i)) = cr(i) 




11 

%0 

i K 


and 

T(v(i)) = a (i) 




ir(i) , 

1 e K 

Now 

T(Tr(k-l)) 

= <J(k-l) 


=X 

T(k) 

= P 



T (Tr(p-l) ) 

= o(p“l) 



T(p) 

= k 



Similarly 


TCk+s+1) = p+t+1 
and T (p+t+1) = k+s+1 
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Thus T is product of two cycles, namely (k,p) and 
(k+s+l,p+t+l). Therefore, by Theorem 3, tt and cr 
are adjacent* 

Above ue have assumed that 1 < k and p+t- < n* 

If k = 1 and p+t < n or k < 1 and p+t = n then there 
is no difficulty involved except that p+t+1 = n+1 
should be taken as 1. However, the case k = 1 and 


p+t = n needs some consideration. We observe that in ! 

this case T(i) i for i = l,p,s+2 only and accordingly 

T is cyclic which implies o is adjacent to tt on 

Tour cf as obtained above 6 Z. or Z.. 

3 

R ule 6 ; Let 

TT = (l,2,3,4^T".7k,T.'7^'3r+l3r+2, . .,r+s^, . . ,r+S 2 j . .,r+Sg, .,r+s^, .,n) 
Let the block elements 1,2,3,4, ...jk, .,«,r be inserted 
in the remaining block of elements 


Cr+X , r+2 , • , • , r+s , • *,r+S2, * , • ,r+s^, ♦ •• , r+ s^ , * • • , n 5 
such that element k is inserted between r+s^ and r+s^+1, 
to give a tour a • Then a is adjacent to tr on 


TP , 
n 

Proof : 

cT (r+1, r+2 , r+o , * » « , r+s^^, 1, r+s^+1, * * , 2 ^+S 2 , 2, r+Sg+l, • * jr+s^, 3, 

r+Sg+l, . .,r+s^^,k, r+s^^+l, . . ,r+s^.,r,r+s^+l, . .,n) 

Then observe that 


2 r+s^^+l 3 r+Sg+1 4 .vr+l r+s^+l 

r+ST+1 1 r+s +1 2 r+So+l..r+s +1 r 
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Clearly is cyclic and the cycle is 

(r+l, r+s^+1 , r , r+S3^]_+l, . . • » 4 , r+s^+l, 3 , r+Sg+l> 2 , r+Sj+1 , 1 ) 
which is of length 2r+l, 

Corol lary 2. ; If more than one element of the block 
(lj2,3, ,, .,k, •. .,r) are inserted between any two elements 
of the remaining block such that the order of occurence 
of the elements (1,2,3 , ,k, jiy) in cf is not 
disturbed i.e, if 3,k £ (1,2, ,,, . .,k, . .,r) and j proceeds 

k then j occurs first in cf than k, then again a is 
cyclic. However, the cycle in this case is not of length 
2r+l. The proof follows obviously by writing cr tT"! as in 
Rule 6. 

Rule 7 : Let rr = (1,2,3, .. •••>n) 

(i) If p = 2s (even) then construct 

a - (1,2,3, .,k-l,l5-l,'bi3^k+5, ... ,k+p-l,k+p,k,k+27fc+47k+i>^.,,,n) 

(ii) If p = 2s+l (odd) then construct 

= (1,2,3, , , . ,k-l,Ic+i, icfsTS^S, • ,775+^2 , k4^, k',fci"2^ k+4 , V.'jlc+p^, . . ,n) 
In •both cases <s is adjacent to tt on TP^^. 

P roof : The details of the proof are omitted here as it 

can be proved on similar lines as the Rule 5 or B\iLe 6. 

However we give its raain points. 

(i) If p = 2s (even) then crir"^ is cyclic and the cycle is of 
length 2s+l. 

(ii) If p = 2s+l(odd) then air“^ is cyclic and the cycle is of 
length 2s+3. 
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As we have seen there can he number of rules 
that can be formulated to generate adjacent tours to 
a given toiar on without exhausting all of them. 

Similarly there can be nuraber of rules that can be 
formulated to generate non— adjacent tours to a given 
tour. Again it is not easy to exhaust all of them 
by formulating such rules. Number of such rules are 
formed in an attempt to separate the set of non-adjacent 
tours to a given tour from the set of tours adjacent to 
it. But always there are some tours left whose adjacency 
or non- adjacency cannot be decided by such systematically 
generated rules. We here do not list all of them but 
give just one rule for non— adjacency which is simple 
in nature and at the same time exhausts a good number 
of to-urs. 

R ule 1 : Let 

- (Tr(l);7r(2_)_ijr(3), . . ^T rCk-l) ,m(k)^ 'Tr(k+l), ««.,Tr(n) ^ be a tour.Let 

be obtained from oc such that ^( 2 ), oO),..., cy(k-l) and 
cf (k+ 1), 0 (n) are some non-identity permutations of 
'ir(2),7 r (3), , . , jir(k-l) and irCk+l), . , .,vCn) respectively. 

Then 3 is non-adjacent to »(• 
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Proof : Construct 

= (ttCI), 0 ( 2 ), 0(3), ...jCf (k~l),Tr(k},B-(k+l), .• ., 7 r(n)) and 

Tg = (Tr(l), v(2)f 'tT(3),,..} Tr(k-l),7r(k),cr (k+1), * • * Cn) ) 

Then observe that 

I x(oC) + I :c(3) =1 :sCT^) + I x(T^) 
which completes the proof. 


6 . 3 * 2 . Applications of Ad.iacency Rules 

In the following, vje describe the applications 
of adjacency rules for the development of heuristic procedures 
for solving the traveling salesman problems, 

Croes [39] , to our knowledge, v/as the first 
to use the idea of permutations for tour to tour improvements 
for the solution of traveling salesman problems* The idea 
of "Inversions” is due to him. Starting with a. initial 
tour TTj an adjacent tour is generated by an inversion, 

(though any such adjacency ideas were not known at that time)v 

'B 

S'' 

If the objective function value for this new tour ("Inverted” 
tour) is less than the previous one, this "inverted" tour is 
retained otherwise a new inversion is tried. Similarly 
there are. some other heuristic methods l22 ] which proceed 
from tour to toin?. It is observelthat in some cases, though 
as a matter of coincidence only, they proceed from a tour 
to an adjacent tour. 
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Reiter and Sherman £22,129] havh described a 
series of fonr similar algorithms that are called ALGOIVCr)# 
In ALGOIV(l), one starts with a random tour, say, 

TT = (i 3 _,i 2 ,i 2 , ,,i^) and finds the best tour amongest 
the follo^-fing tours 

^ ^ , i^ , • * • , iL ) 

<j ~ (i_, 1 j i j • » • , i ) 

2 2 3 1 n 

« 

* 

Vs " Vl’h’V 

(i.e, find the best position to insert i^ in the sequence 

(igjigji^j ••• ,ij^))« It can be observed by Rule 6 or by Rule 1 

of cyclic changes that all such tours are adjacent tours* 

However, one is testing here only (n-2) neighbouring 

n~l 

vertices while cyclic changes generate 2 (r-l ) (n+l-r ) 

r=2 

neighbouring vertices, a number much larger than (n-2). 

If tour TT - (i^jigjig, ...,i^) then ALG0tTr(2) 

seeks a best possible insertion of 'the sequence 

(i^,i 5 ic.j*..ji ). Although the generated! tours are not 
in accordance with Rule 6 (since in the insertion of 
i„ may proceed i-) but still they are adjacent as 
air”^(i) i for at most five values of i. 

Similarly ALG0IV(3) seeks a best possible insertion 
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of ijQ the sequence (assuming 

the current tour as (i-, ,i_,i_,i. ,i£.,i^, *.,,i )), But 
ALG0I?(3) always does not generate adjacent tours e,g* 
if IT = (T*2 3 4 5 6 7 8) and the tour generated 

is O' = (4 2 5 1 3 6 7 8) then <3t is non-adjacent 

to IT, Similarly 1LG0IV(4) generates some adjacent and 
some nonr-adjacent vertices.. 

Reiter and Sherman has coupled all the four algorithms 
together. When ALGOIV(l) fails to seek an improvement, the 
control is transferred to ALG0IV(2). After the improvement, 
control again returns to ALGOlV(l). If ALGOIV(l) and 
ALG0IV(2) both fail to give an improvement, the control 
is transferred to ALG0I¥(3)* After the improvement, 
control again returns to ALGOIV(l). Thus when ALG0IV(3) 
fails to give an Improvement, control is referred to 
ALG0IV(4)* The tour is declared optimal when ALG0IV(4) 
fails to give an improvement. While ALG0IV(3) and 
ALG0I7(4) do generate some non~adJacent vertices also, 
ALG0I¥(1) and ALG0IV(2) to which control is referred 
mainly generate adjacent vertices only. However, there 
was no such idea of adjacency of tours even at the time 
of Reiter and Sherman [ 129 ] , 

Thus it adds to our belief that by using adjacency 
rules or the combinations of a few of them efficient heuristic 
procedures can be developed for solving traveling salesman 
problems. 
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A graph is said to he Hamiltonian if it contains 
■ a circuit x^hich includes all nodes of the graph* Let 
G(TP^) be the graph of defined similar to graph 
G(AP^) in Chapter 5, 

Theorem 6 : G(TP^) is hamiltonian* 

Proof ; We knov/ that the interchange of any two consecutive 
elements of a tour r results in an adjacent tour 
(Rule 1), Using .this rule and following the proof given 
hy Balinski [12] the proof can be completed* 

^•4 Further Scone of Research 

We here list some problems for further 
research. A positive answer to these problems will 
contribute a lot to the traveling salesman polytope 
and might help in developing some algorithm for the 
traveling salesman problem also* 

(1) To find necessary and sufficient conditions for 
two tours to be adjacent (or non- adjacent) on TP^. 

(2) To find the number of tours adjacent (on 

to given tour, Hox^ever, fair evidence is there that such a 
nuiaber is going to be a very large number* 

(3) For a given tour v and c such that t and 
are non- adjacent on TP , Padberg and Rao [l2l] have 
given a method to construct a tour T such that w-i and 
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T ^ are cyclics and accordingly T is adjacent to both 
hr md. a on AP^ and thus also on TP • However, a tour T 

m -1 

can be adjacent to both it and u without being Tir 
fUid as cyclics. To exhaust the set of all those 

tours T such that T is adjacent to both tt and cr 
on TP^ is another associated problem, (By taking w, c 
find T as assignments one c£xn think of similar problem 
in case of AP^^). 

(4) If w and a are two adjacent tours on then 

we know = v where v is cyclic i.e* o = vtt can be 

generated from tt by the knowledge of v* Thus the problem 
is to find out all those possible v’s (cyclics) such 
that generated cj»s are tours* Apparently some conditions 
on V are needed to exhaust all those v’s which generate 
tours adjacent to tt and AP^. Rao [l26l has given one 
such V to construct ^ adjacent to Tr( on AP^^). 

One condition on v is that v must have odd number 

1 

of cleuionts * This is clear from the fact that v = cfT^ 
is an even permutation, 

(5) Let TT be the optimal assignment which is not a 
tour. Let B(ir) be the set of all those tours which are 
adjacent to ir on AP^, Then the problem is to determine 
the set BCtt), A few members of BCtt) are determined belo^ 

Let TT be the product of k disjoint cycles (cycles 
of length one are also included), We know, as shown by 
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Baliiiski tlsJ , also sea section 5.4, that w can "be 
written as product of two permutations one of which 
is a tour and the other is cycle of length k* 

Thus TT = uv where u is a tour and v is cycle of 
length k. Clearly u~^ = v is cyclic which implies u 
is adjacent to tt on AP^. If r^^jr^, . . . ,r^ are the 
lengths of k cycles of ir then possibles choices of 
u( or v) are 

(r j^r «. C h"" 1 ) 1 

which gives minimum number of elements of BCtt). The 
following result helps to some extent to know about 
the nature of members of BCtt). 

Re suit i If c 8 BCtt) then a 8 Z (tt) for some p > k. 

Jr mm 

Because, otherwise o it"! is cyclic of length 
p < k which implies tt and o have at least one cycle 
common with each other. This contradicts that cr 
is a tour. Hence p > k. 
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